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ABSTRACT 

A  digital  data  acquisition  system  was  implemented  for  the  study  and  collection 
of  scattered  waveforms  from  various  simple  objects.  Reflected  waveforms  were 
extracted  from  interference  with  other  signals  by  a  simple  subtraction  technique. 
Backscattered  and  bistatic  reflections  were  observed  and  recorded.  Various  effects 
on  scattered  waveforms,  including  target  composition,  pulse  length,  and  frequency, 
were  experimentally  observed.  Bistatic  scattering  results  from  a  stryofoam  sphere 
at  a  ka  value  of  10  were  in  agreement  with  theoretical  predictions  for  a  rigid  sphere. 
Bistatically  measured  data  at  70  kHz  from  a  styrofoam  cylinder  of  length  20  cm  and 
diameter  2.5  cm  followed  a  trend  for  bistatic  angles  up  to  150  degrees  but  became 
amorphous  as  the  bistatic  angle  approached  180  degrees. 
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I.  TNTRODUCTION 

Any  object  with  a  characteristic  impedance  different  from  the  surrounding 
fluid  will  affect  an  acoustic  wave  impinging  on  it.  The  wave  will  be  transmitted 
through  the  object  and  scattered  from  it.  Flexural  waves,  compressional  waves,  and 
creeping  waves  are  examples  of  some  of  the  possible  scattering  and  reflection 
mechanisms  [Ref.  l:p.  301].  The  "characteristic  impedance"  equals  the  density  of  a 
material  multiplied  by  the  speed  of  sound  in  the  material.  It  is  one  measure  used  to 
calculate  reflection  factors  for  pressure  waves.  Other  properties  that  affect  the 
scattering  properties  of  an  object  include  frequency,  shape  of  the  object,  and  pulse 
length. 

Ever  since  the  development  of  the  submarine  as  an  important  force  in  naval 
warfare,  much  effort  has  been  directed  to  improving  their  detection.  The  same 
holds  true  for  the  submarine  detecting  of  surface  ships  or  other  submarines.  The 
knowledge  required  for  detection  include  knowing  the  source  level  of  the  target, 
propagation  loss,  noise  level ,  and  lastly  target  strength,  the  subject  of  this  paper. 
These  parameters  are  used  in  the  well-known  sonar  equations  [Ref.  2:p.  411]. 
Target  strength,  TS,  as  used  in  the  sonar  equations  is  defined  as: 

TS  =  101ogi^  (1) 

where  If  is  the  reflected  intensity  one  meter  from  the  acoustic  center  of  the  target 
and  lo  is  the  incident  intensity  at  the  target  without  the  presence  of  the  target.  This 
equation  is  general,  giving  no  stipulation  as  to  the  direction  of  the  reflected 
intensity.  Typically,  signals  in  active  sonar  are  sent  and  recevied  by  the  same 
platform,  thus  the  most  common  geometry  and  ,therefore,  the  most  studied,  is  the 


backscattered  monostatic  case  where  the  reflected  intensity  is  measured  in  the 
direction  back  toward  the  source.  Measurement  along  any  other  direction 
constitutes  a  bistatic  geometry.  Relatively  little  effort  has  been  given  to  the  study  of 
bistatic  target  strength  until  recently  when  the  bistatic  geometry  was  recognized  as  a 
potential  method  for  improving  underwater  detection  of  submerged  targets. 
Figures  1  and  2  illustrate  the  monostatic  and  bistatic  geometries. 
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Figure  1.  Backscattered  Monostatic  Geometry 


Source 


Target 


Receiver 


Figure  2.  Bistatic  Geometry 


Some  work  has  been  done  in  the  radar  community  on  estimating  the  bistatic 
radar  cross-section  of  objects.  Crispin  and  Siegel  [Ref.  3:p.  156]  present  a  theorem 
for  the  approximation  of  radar  bistatic  cross-section: 

In  the  limit  of  vanishing  wavelength,  the  bistatic  cross-section  for 
transmitter  direction  tc  and  receiver  n  is  equal  to  the  monostatic  cross-section 
for  the  transmitter-receiver  direction  k  -i-  n  with  k  ^  n  for  bodies  that  are 
sufficiendy  smooth. 

(Figure  3  illustrates  the  geometry  as  stated  in  the  theorem.) 
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Figure  3.    Radar  Bistatic  Cross-Section  Theorem  Geometry 

Their  bistatic  theorem  has  its  limitations  wHen  applied  to  complex  shapes;  The 
bistatic  cross-section  can  be  greater  or  less  than  the  monostatic  case  due  to 
structurally-related  interferences.  Another  limitation  is  the  size  of  the  target  with 
respect  to  the  incident  wavelength.  Spheres  were  analyzed  using  various  ka  values 


where  where  a  is  the  radius  of  the  sphere  and  k  is  the  wave  number.  It  was  found 
that  the  theorem  held  more  closely  for  higher  ka  values,  but  in  all  cases  the  theorem 
sharply  deviated  from  experiment  as  the  bistatic  angle  approached  180  degrees.  As 
an  illustrative  example,  they  presented  data  on  a  jet  aircraft  and  a  ballistic-missile 
shape  [Ref.  3:p.  174-177].  Data  were  in  fairly  good  agreement  with  their  theorem 
as  long  as  the  bistatic  angle  did  not  get  too  large. 

In  transition  from  the  radar  community  to  underwater  acoustics,  Urick  [Ref. 
l:p.  302]  implied  that  the  bistatic  radar  cross-section  theorem  might  apply  in  a 
underwater  condition,  in  which  case  the  bistatic  target  strength  was  equal  to  the 
monostatic  target  strength  taken  at  the  bisector  of  the  bistatic  angle.  Urick 
emphasized  there  was  a  lack  of  data  to  prove  or  disprove  such  a  theory.  It  would  be 
very  convenient  if  such  a  theorem  existed  since  there  is  plentiful  target  strength  data 
for  monostatic  geometries.  Whether  or  not  the  high  frequency  limitations  of  such  a 
theorem  would  preclude  its  tactical  use  is  another  question. 

The  interest  in  bistatic  target  strength  has  been  generated  by  the  possibility  of 
increasing  detection  of  submerged  and  surface  units  using  this  technique.  For 
example,  one  can  easily  visualize  having  a  surface  ship  using  its  active  sonar  while  a 
submerged  submarine,  using  passive  sonar,  listens  for  reflections  off  enemy  forces. 

There  are  several  advantages  and  disadvantages  to  bistatic  employment.  Several 
ideas  now  discussed  were  taken  from  Harvey  [Ref.  4]  who  gave  a  thorough 
discussion  on  using  bistatic  employment  as  a  tactical  tool  in  naval  warfare.  The 
major  advantages  include  the  ability  to  maintain  a  track  of  the  enemy  force  using 
passive  sonar  and  the  increased  amount  of  data  available  to  detect  and  classify 
targets.  All  units  besides  the  active  unit  can  be  passive,  thereby  hampering  own- 
ship  detection.  Detecting  bistatic  reflections,  as  well  as  backscattered  reflections, 


from  a  target  increases  the  data  available  for  detection  and  classification  of  the 
target.  The  greatest  difficulty  using  a  bistatic  geometry  is  the  high  positional 
accuracy  required  between  the  pinger  and  the  listener.  This  requires  prior  planning 
for  use  with  surface  ships  and  submarines  since  they  do  not  normally  have 
continuous  open  communication  with  each  other.  Another  difficulty  is  the 
distinguishing  of  backscattered  reflections  from  bistatic  reflections.  It  may  turn  out 
that  the  advantages  do  not  outweigh  the  disadvantages  and  the  whole  idea  will  prove 
infeasible. 


II.  TARCwET  STRENGTH 

A.  INTRODUCTION 

This  chapter  presents  some  definitions  and  background  knowledge  essential  in 
preparation  for  experimentation.  Some  concepts  are  rather  simple  but  are  included 
for  completeness. 

B.  DEFINITION 

The  definition  of  target  strength  has  already  been  given,  equation  (1). 
However,  there  are  several  aspects  of  the  definition  that  warrant  clarification. 
First,  all  measurements  are  taken  in  the  far  field  of  the  target.  The  far  field  is 
where  the  pressure  reflected  from  an  object  exhibits  a  1/r  relationship  [Ref.  2:p. 
178].  For  example,  pressure  measurements  taken  one  meter  from  the  hull  of  a  large 
ship  would  certainly  be  in  the  near  field,  so  intensity  is  measured  in  the  far  field  at 
ranges  of  hundreds  of  meters  and  then  extrapolated  back  to  one  meter  for  use  in  the 
target  strength  calculations. 

The  next  question  concerns  the  defintion  of  the  acoustic  center.  Very  rarely  are 
we  able  to  know  the  location  of  a  target's  acoustic  center.  The  acoustic  center  is  the 
point  from  which  sound  appears  to  radiate.  It  may  be  anywhere  around  the  target, 
and  may  lie  outside  the  shape  itself.  It  may  also  change  with  the  aspect  of  the  target 
relative  to  the  receiver.  Exact  location  of  the  target's  acoustic  center  is  ignored  in 
field  measurements,  and  the  geometrical  center  is  used.  When  measuring  pressure 
at  hundreds  of  meters,  a  few  meters  error  in  determination  accounted  for  by  the 
acoustic  center  will  not  produce  a  large  discrepancy. 


C.    SPHERICAL  WAVES 

For  a  distant  source,  to  a  good  approximation,  the  wave  will  appear  locally  as  a 
plane  wave  as  it  encounters  the  target.  Consider  an  omnidirectional  point  source: 
Spherically  diverging  waves  are  emitted  and  the  pressure  amplitude  and  phase  is 
dependent  only  on  the  distance  r  from  the  acoustic  center  of  the  source.  An 
outwardly  travelling  spherical  wave  can  be  represented  by 


o/  X     A     J  (wt  -  kr) 

P(r)=— e  (2) 

r 


A  =  pressure  amplitude  at  one  meter 
r  =  distance  from  point  source 
w  =  angular  frequency  of  signal 
k  =  wave  propagation  number 


When  the  radius  of  curvature  of  the  wave  front  is  large  compared  to  the  biggest 
dimension  of  the  target,  the  spherical  wave  can  be  treated  as  a  plane  wave. 

Plane  waves  are  assurred  when  the  dimensions  of  a  target  are  small  compared  to 
the  source-target  distance.  One  of  the  goals  of  small-scale  experiments  is  the 
application  of  the  results  to  oceanic  situations,  therefore  there  is  a  concern  about  the 
approximation  of  plane  waves  ensonifying  a  target  in  a  small  laboratory  tank. 

The  geometrical  constraints  of  the  water  tank  limits  target-source-receiver 
distances.  With  small  distances  between  source  and  receiver  and  big  targets,  a  plane 
wave  may  not  be  a  good  approximation.  Another  factor  that  influences  the  plane 
wave  approximation  is  the  directionality  of  the  source.  It  is  difficult  to  achieve 
omnidirectional  behavior  from  a  finite  source  and  it  is  impossible  if  the  frequency 
is  high.    Beam  patterns  are  formed  with  energy  emitted  in  preferred  directions.  If 


the  target  subtends  an  angle  greater  than  the  beamwidth,  the  pressure  amplitude  will 
not  be  the  same  for  all  points  on  the  target. 

Plane  waves  also  have  properties  that  make  target  strength  calculations  simple 
to  handle.  The  intensity  can  be  equated  to  the  root-mean-square  pressure  amplitude 
squared  divided  by  the  characteristic  impedance  [Ref.  2:p.  111].  Since  hydrophones 
measure  pressure  via  a  sensitivity  factor,  the  output  voltage  they  produce  can  be 
directly  related  to  intensity. 

D.    WAVEFORM  EFFECTS 

There  are  several  factors  to  consider  in  the  measurement  and  collection  of 
reflected  waveforms  from  objects.  Some  are  listed  below  and  each  will  cause  the 
reflected  waveform  to  differ  from  the  incident  waveform. 

-  target  composition 

-  target  shape 

-  pulse  length 

-  resonance  effects 

For  the  reflected  waveform  to  be  identical  of  the  incident  waveform,  a  planar 
rigid  boundary  must  be  encountered.  For  more  complicated  targets,  some  of  the 
incident  energy  will  be  transmitted  into  the  body  and  eventually  reradiated  and 
different  areas  of  the  target  may  reflect  in  such  a  way  to  cause  mutual  interference. 
Many  objects  have  resonances,  at  which  the  target  will  store  some  of  the  incident 
energy  and  reradiate  it  at  a  later  time  causing  interference  effects.  Theoretical  and 
experimental  results  show  that  for  short  pulse  lengths,  reflections  are  nearly 
independent  of  target  composition  and  are  almost  identical  to  the  incident  pulse 
[Ref.  5,6,7]  while  for  longer  pulses  the  reflected  waveform  changes  drastically 
depending  on  the  composition  and  shape  of  the  target. 
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E.    EXPERIMENTAL  MEASUREMENT 

While  conducting  the  literature  search,  it  was  interesting  to  note  how  the 
pressure  amplitude  of  the  reflected  wave  was  measured.  As  stated  before, 
hydrophones  give  a  voltage  that  is  representative  of  pressure  and  this  can  be  equated 
to  intensity.  The  problem  lies  in  the  fact  that  the  reflected  waveform,  in  most  cases, 
is  not  identical  to  that  of  the  incident  signal.  The  experimenter  has  to  specify  what  is 
measured  and  what  kind  of  pulse  is  being  used.  Response  can  be  catagorized  as 
transient  or  steady-state.  Steady-state  implies  that  the  whole  target  has  been 
ensonified  by  the  incident  pulse  allowing  for  all  portions  of  the  target  to  contribute 
to  the  reflection  and  for  any  resonances  to  become  fully  excited.  In  this  case,  the 
return  will  be  a  combination  of  all  possible  reflection  paths  and  modes.  A  transient 
response  simply  is  not  in  the  steady-state.  To  give  an  idea  of  measurement  used,  the 
following  are  quoted  from  two  published  articles: 

...five  cycles  of  the  sine  wave  located  at  a  position  in  the  pulse  judged  to 
be  a  true  representation  of  the  steady  state.  The  average  of  the  ten  peak  values 
so  obtained  gave  the  peak  pressure  amplitude  [Ref.  7] 

After  amplification,  the  amplitude  of  the  back- scattered  signal  is 
measured  with  a  peak  detector  [Ref.  8]. 

The  main  point  is  that  for  reflected  waveforms,  not  similiar  in  shape  to  the 
incident  waveform,  some  subjective  determination  is  required.  Estimation  of  the 
steady-state  must  be  made  or  some  definite  parameter  has  to  be  selected  such  as  peak 
amplitude  or  total  energy.  In  any  case,  an  unique  way  to  convert  the  amplitude 
measurements  to  an  intensity  level  is  not  achievable.  Add  to  this  the  fact  that  the 
reflected  waveform  will  be  different  for  different  target  orientations  and 
frequencies,  and  the  problem  of  translating  measurements  into  an  intensity  level  to 
calculate  a  target  strength  becomes  complicated. 


There  presently  exists  a  Target  Strength  Standardization  Program  [Ref.  9] 
which  was  assigned  the  task  of  finding  some  standard  measure  of  target  strength 
independent  of  signal  type.  The  general  consensus  from  the  program  was  to  use  the 
"Average  Integrated  Target  Strength"  which  is  the  average  time- integrated  energy. 
Furthermore,  an  Echo  Metrics  Subcommittee  [Ref.  9]  was  formed  to  find  signal- 
processor  dependent  metrics  for  the  measurement  of  target  strength.  Theoretically, 
the  incident  energy  encountered  by  the  object  would  be  reflected  or  absorbed  and 
eventually  reradiated  at  some  later  time  providing  no  energy  was  stored  or  lost  in 
the  object.  Therefore,  it  would  seem  that  the  average  of  the  integrated  energy 
would  be  an  independent  metric  for  all  types  of  signals  providing  the  integration 
time  is  extended  to  include  all  the  return  signal.  Some  standardization  of  target 
strength  is  being  investigated  so  that  data  from  independent  experiments  can  be 
compared  and  integrated. 
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III.   DATA  ACOTJTSTTTON  SYSTEM 

A.  BACKGROUND  AND  CONCEPT 

The  experiments  were  conducted  in  a  water  tank  located  in  Room  019  in 
Spanagel  Hall  at  the  Naval  Postgraduate  School.  The  tank  was  approximately  1.5  m 
wide,  8  m  long,  and  2  m  deep  and  lined  with  anechoic  rubber  material.  For  the 
targets  used,  wall  reflections  were  of  the  same  magnitude  as  the  target  reflections. 
Usually,  it  was  easier  to  use  pulsed  signals  and  time  discriminate  for  the  target 
reflection.  However,  the  small  dimensions  of  the  tank  did  not  provide  the 
environment  to  do  this  effectively  for  all  geometries.  To  overcome  this  difficulty,  a 
subtraction  technique  was  used  employing  a  digital  storage  device.  The  basic 
concept  was  to  subtract  the  signal  received  with  the  target  immersed  in  the  tank 
from  the  signal  received  with  the  target  absent.  The  procedure  was  very  simple  and 
allowed  for  any  target-source-receiver  geometry. 

B.  DEVELOPMENT 

A  considerable  amount  of  time  went  into  the  selection  of  the  digital  sampling 
and  storage  device,  the  source,  and  the  receiver.  The  first  question  that  had  to  be 
addressed  was  the  range  of  frequencies  to  be  used.  The  frequency  range  of  100-200 
kHz  was  initially  selected  because  it  gave  wavelength-to-target  dimension  ratios  that 
match  oceanic  scenarios.  The  critical  specifications  for  the  digital  storage  device 
are  the  sampling  rate,  dynamic  range  and  jitter.  Dynamic  range  was  not  a  problem 
since  most  digital  oscilloscopes  have  the  same  voltage  range  as  a  typical  non-digital 
oscilloscope.  The  sampling  rate  must  be  fast  enough  to  provide  a  good 
representation  of  the  signal,  and  assuming  a  sinusoidal  signal  about  ten  samples  per 
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cycle  would  give  a  good  trace.  For  example,  for  a  100  kHz  signal  the  sampling  rate 
would  have  to  be  1  MHz.  Another  major  concern  was  the  amount  of  jitter  in  the 
sampling  device.  Jitter  describes  the  fluctuation  in  the  triggering  times  relative  to 
the  sampled  signal.  For  a  system  using  a  subtraction  process,  it  is  perhaps  the  most 
critical  specification  because  there  has  to  exist  a  common  time  reference  between 
collected  waveforms  to  get  accurate  results. 

The  selection  of  the  digital  oscilloscope  was  based  on  the  availability  of 
equipment.  Two  digital  oscilloscopes  were  evaluated.  The  Nicolet  3091  was 
analyzed  first.  Its  fastest  sampling  rate  was  1  MHz  and  had  a  storage  capability  of 
4000  points.  The  Nicolet  3091  had  a  RS-232  serial  interface  and  was  easy  to  use, 
however  it  was  not  acceptable  due  to  excessive  jitter  at  the  frequencies  of  interest. 
The  jitter  was  easily  seen  by  applying  a  pulsed  signal  and  observing  the  traces  on  the 
scope.  The  digital  oscilloscope  finally  used  was  the  Tektronix  2430  (Tek  2430).  It 
is  a  very  capable  oscilloscope  and  has  many  nice  features  making  it  suitable  for  data 
acquisition.  The  most  important  feature  that  it  had  was  an  averaging  mode 
selectable  from  2  to  256  waveforms.  This  helped  eliminate  random  noise  and 
minimized  any  fluctuations  the  electronic  equipment  possessed.  The  amount  of 
jitter  was  dependent  of  the  time  scale  selected.  The  fastest  sampling  rate  of  the  Tek 
2430  was  100  MHz  which  could  be  used  to  sample  a  10  MHz  signal  using  ten  points 
per  cycle.  The  jitter  specification  given  in  the  operating  manual,  for  the  time- 
delayed  mode,  was  .08  x  (time  scale).  The  jitter  was  observed  and  was  significant 
for  the  time  scale  desired,  however  the  saving  factor  was  the  time  averaging  mode 
which  ultimately  averaged  out  the  jitter  and  provided  a  sufficiently  constant  time 
reference.  It  was  very  easy  to  observe  this  on  the  oscilloscope  because  of  a 
triggering  marker  on  the  signal.   Jitter  becomes  less  of  a  problem  for  any  digital 


12 


oscilloscope  when  the  number  of  samples  per  cycle  of  signal  increases.  However, 
the  problem  then  arises  from  the  limited  amount  of  storage  in  the  oscilloscope.  A 
smaller  number  of  cycles  can  be  stored  when  using  a  faster  sampling  rate.  The  Tex 
2430  has  a  storage  capabiltiy  of  1024  points,  much  less  that  the  Nicolet  3091; 
however  it  also  has  a  time-delay  feature  which  allowed  for  selective  data  storage 
(The  Nicolet  3091  did  not  have  any  time-delay  feature).  Other  features  of  the  Tek 
2430  include  cursor  control  for  voltage  and  time  measurements,  simple  output  to  a 
Hewlett  Packard  Think- Jet  printer  or  pen  plotter,  and  a  GPIB  interface.  It  is  a  very 
versatile  piece  of  equipment  and  was  finally  selected  for  use  in  the  data  acquisition 
system. 

The  selection  of  the  source  and  receiving  hydrophone  was  also  limited  to  in- 
house  equipment.  For  the  frequencies  desired,  the  USRD  type  F41  transducer  was 
the  appropriate  choice  with  a  nominal  frequency  range  of  15-150  kHz  [Ref.  10:p.  5]. 
It  is  a  piston  type  piezeoelectric  transducer  with  a  face  diameter  of  10  cm.  The 
choice  of  the  receiver  was  a  Celesco  LC-10  hydrophone.  It  has  a  wide  nominal 
frequency  band,  is  small  and  light  making  it  easy  to  move  around,  and  has  good 
sensitivity. 

C.    SYSTEM  SET-UP 

1.    Electronic   Equipment 

The  F41  source  had  to  be  driven  near  200  V  rms  for  target  reflections  to  be 
seen.  An  Instruments,  Inc.  Model  L2  power  amplifier  was  used  to  drive  the  source 
and  its  selection  was  again  based  on  what  was  available  in  the  laboratory.  It  allowed 
for  impedance  matching  and  has  a  voltage-divided  10  mVfV  port  to  monitor  output 
voltage.  The  output  of  the  receiver  was  passed  through  an  Ithaco  Model  1201  Low 
Noise  preamplifier  which  has  bandwidth  setting  capability.    This  was  very 
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beneficial  because  there  was  a  lot  of  60  Hz  noise  which  had  to  be  eliminated  by 
bandwidth  filtering.  The  signal  out  of  the  preamplifier  was  sent  to  a  Kikusui 
COS5060A  oscilloscope  and  the  Tektronix  2430  digital  oscilloscope.  Having  a  non- 
digital  oscilloscope  made  observing  and  measurement  of  the  waveforms  easier  to 
perform  on  the  Tek  2430.  A  Hewlett  Packard  3314A  function  generator  was  used 
to  provide  the  source  signal.  It  also  provided  the  triggering  signal  which  was  sent  to 
the  external  channels  of  the  oscilloscopes.  In  the  pulse  mode,  the  HP  3314A  sends 
out  a  3  V  transistor-transistor  logic  signal  via  its  trigger  port.  The  pulse  repetition 
rate  of  10  ms  was  sufficient  to  allow  reflections  to  diminish  before  the  next  pulse.  A 
Hewlett  Packard  9000  Series  300  desk  top  computer  was  used  as  the  data  collection 
and  processing  unit.  It  was  connected  to  the  Tek  2430  via  the  HPIB  bus.  The 
transfer  of  waveforms  stored  on  the  Tek  2430  was  controlled  by  software 
programs.  Data  was  ASCII  encoded  from  the  Tek  2430,  transferred,  decoded  into 
actual  voltages,  and  then  stored  onto  disks  in  binary  data  form.  A  schematic 
diagram  of  the  electronic  system  is  shown  in  Figure  4. 
2.    Physical  Set-uo 

To  measure  all  possible  bistatic  angles,  the  scattering  pattern  would  have  to 
obtained  in  three  dimensions.  This  was  not  done  in  this  work.  Bistatic  data  were 
limited  to  the  horizontal  plane  determined  by  the  source,  receiver,  and  target,  with 
the  axis  of  the  target  in  the  same  plane. 

The  source,  the  F41  transducer,  was  mounted  on  a  3/4-inch  aluminum  rod 
attached  to  a  rectangular  support  plate.  The  rod  was  fixed  to  the  plate  by  means  of  a 
friction  fitting  secured  by  a  hose  clamp.  Target,  source,  and  receiver  were 
submerged  one  meter  below  the  surface  of  the  water  and  the  rods  which  supported 
them  were  marked  for  this  depth.   The  source  with  its  support  rested  on  another 
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Figure  4.  Data  Acquisition  System 


15 


support  built  from  aluminum  L-beams  that  spanned  the  width  of  the  tank.  Lead 
bricks  were  used  to  hold  the  source  once  it  was  set  in  place  and  a  level  was  used  to 
ensure  the  source  face  was  vertical.  The  receiver  was  mounted  onto  a  sliding  track 
specially  built  for  the  tanks.  Adjusting  screws  were  added  for  ease  of  positioning 
the  receiver  vertically.  The  initial  intent  was  to  design  a  rotating  boom  to  hold  the 
receiver  and  support  the  target  from  its  center.  A  device  was  constructed,  but  too 
much  receiver  movement  occurred  between  the  loading  and  unloading  of  targets. 
Though  the  boom  was  never  used,  the  device  was  still  used  for  the  suspension  of 
targets  and  had  a  mechanical  angle  indicator  on  its  top  for  angle  measurement. 
Spherical  targets  were  suspended  by  a  thin  cable  and  the  cylindrical  target  by  a  3/8- 
inch  rod. 

The  direction  of  the  source  was  set  by  aligning  the  source,  receiver,  and 
target  along  a  straight  line  with  the  receiver  between  the  source  and  target.  Source, 
receiver,  and  target  depth  were  one  meter.  Vertical  alignment  was  ensured  by 
keeping  the  source  and  receiver  vertical  by  use  of  a  level.  The  frequency  was  set  to 
200  kHz  to  establish  a  narrow  major  beam  and  the  source  was  rotated  until 
maximum  signal  was  seen  by  the  receiver.  This  ensured  the  target  was  being  totally 
ensonified  by  the  major  lobe  of  the  source. 

Once  the  source  direction  was  set,  the  receiver  was  moved  to  the  desired 
position.  The  angle  between  the  source,  target,  and  receiver  was  determined  by 
setting  the  receiver's  track  perpendicular  to  the  source-target  axis  and  then  using 
trigonometry  relationships.  After  the  desired  geometry  was  set,  care  was  taken  to 
ensure  the  target,  source,  and  receiver  were  stationary.  It  was  critical  that  all 
elements  be  stationary  when  waveforms  were  collected. 


16 


3.     Data   Collection 

The  method  of  data  collection  and  storage  followed  the  sequence: 

-  set  up  desired  source-receiver-target  geometry 

-  warmup  all  electronic  equipment  >1  hr 

-  ensure  connections  as  illustrated  in  Figure  5 

-  load  "Tex"  onto  computer 

-  set  up  the  Textronix  2430  oscilloscope 

-  B  trigger  mode  to  "runs  after" 

-  A  trigger 

-  mode:  normal 

-  coupling:  AC 

-  source:  extemal  1 

-  slope:  positive 

-  level:  800  mV 

-  trigger  position:  1/8 

-  data  encoding:  ASCII 

-  GPIB  address:  4 

-  vertical  mode:  channel  2 

-  select  time  delay  to  obtain  desired  time  window 

-  set  acquisition  mode  to  256  averaging 

-  adjust  channel  2  voltage  scale/Ithaco  gain  appropriately 

for  the  signal 

-  for  each  new  signal,  press  the  acquire  button 

-  wait  sufficient  time  for  averaging  to  take  place 

-  press  the  save  button 

-  run  the  "Tex"  program 

The  "Tex"  program  controlled  the  transfer  and  storage  of  data  from  the 
digital  oscilloscope.  Inputs  from  the  user  were  file  names  and  file  descriptions. 
The  transfer  process,  including  writing  onto  the  disk,  took  approximately  10  s  per 
waveform.  A  listing  of  the  program  is  in  Appendix  A. 
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The  computer  program  "Plot"  was  used  to  evaluate  the  stored  data.  It 
took  two  files  and  subtracted  them,  and  then  plotted  the  results  on  the  CRT 
terminal.  Hard  copy  output  could  be  obtained  by  porting  to  a  pen  plotter  or  Think- 
Jet  printer.  Input  from  the  user  required  file  names  and  where  printout  was  to 
occur.  The  plotting  routine  was  written  so  that  the  plot's  parameters  were  easily 
changeable.  A  hsting  is  enclosed  in  Appendix  B. 

The  method  of  data  collection  was  a  slow  and  somewhat  tedious  process 
since  all  movement  of  target  and  receiver  ware  manually  controlled.    Data 
readings  were  also  recorded  manually. 
4.    Validation 

The  question  of  whether  the  data  acquisition  system  could  correctly 
extract  a  target  reflection  immersed  within  other  signals  had  to  be  answered.  The 
very  simple  procedure  of  subtracting  background  noise  works  only  if  there  is  a 
common  time  reference  for  "target  in"  and  "target  out"  waveforms.  For  ease  of 
discussion,  "background  noise"  is  used  for  any  signal,  random  or  coherent,  that  is 
not  the  desired  signal  and  corresponds  to  the  "target  out"  condition. 

The  procedure  to  validate  the  system  began  with  obtaining  a  geometry 
which  allowed  a  direct  pulse  and  reflected  pulse  to  be  completely  separated  by 
time  discrimation  with  no  interfering  wall  reflections.  Figure  5  is  such  an  example. 
The  larger  pulse  was  the  direct  path  pulse  and  the  second  pulse  is  the  pulse 
reflected  from  a  styrofoam  sphere.  The  frequency  was  66.5  kHz  and  the  pulse 
length  was  20  cycles.  The  reflected  pulse  is  shown  enlarged  in  Figure  6.  The 
peak-to-peak  voltage  of  the  reflected  pulse  was  measured  to  be  68  +  2  mV  using 
the  cursor  controls  of  the  Tek  2430. 
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The  pulse  length  was  then  increased  to  120  cycles  so  the  direct  pulse  would 
be  interfering  with  the  reflected  pulse.  No  other  changes  were  made.  The 
interference  of  the  two  signals  is  shown  in  Figure  7  and  was  stored  in  file  A.  Next, 
the  sphere  was  removed  and  the  background  noise  sampled  and  recorded  in  file  B. 
The  background  noise  in  this  case  was  mostly  the  direct  pulse.  The  last  step  was  to 
subtract  file  B  from  file  A  to  get  the  reflected  signal  from  the  sphere.  The  result  is 
shown  in  figure  8.  The  waveform  begins  at  the  same  time  in  comparison  with 
Figure  6  and  the  only  difference  is  that  the  pulse  length  was  longer  as  expected  and 
the  curve  slightly  more  choppy  in  form.  Measured  voltage  was  73  +  4  mV.  The 
larger  uncertainty  arose  due  to  the  combination  of  small  jitter  and  amplitude 
fluctuations  of  the  HP  3314A.  The  uncertainty  was  estimated  beforehand  by  taking 
identical  signals  directly  from  the  function  generator  to  the  Tek  2430,  recording 
them  in  separate  files,  and  subtracting  them  and  observing  the  results.  It  was  found 
that  peak-to-peak  voltage  uncertainty  was  approximately  four  digitizing  levels 
regardless  of  the  voltage  scale  used  (a  digitizing  level  equals  the  voltage  per  division 
divided  by  50  for  the  Tek  2430).  Given  that  this  error  would  always  be  present, 
there  was  good  agreement  between  Figure  6  and  Figure  8,  and  the  method  of 
obtaining  data  was  considered  validated. 

As  with  any  signal  processing  system,  a  high  signal-to-noise  ratio  is 
desired.  For  any  digital  subtraction  routine,  the  desired  signal  has  to  be  large 
enough  relative  to  the  noise  so  it  will  not  disappear  in  the  digitizing  level  of  the 
sampling  device.  For  example,  if  noise  plus  the  signal  were  10  volts  using  512 
digitizing  levels  or  bins,  each  bin  represents  19.5  mV.  If  the  signal  were  small,  i.e., 
1  mV,  the  signal  level  would  be  much  smaller  that  the  digitizing  level  and  would  be 
lost  in  the  subtraction  routine. 
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Figure  5.  Direct  and  Reflected  Pulse 
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Figure  6.  Reflected  Pulse 
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Figure  7.  Interference  Between  Direct  and  Reflected  Pulse 
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Figure  8.  Reflected  Pulse,  N=120 
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Keeping  the  source  and  receiver  positions  fixed  relative  to  each  other  was 
critical.  Small  movement  between  the  two  greatly  affected  the  phase  of  the  received 
signal  and  had  a  large  effect  on  the  results.  This  effect  was  verified  in  preliminary 
work.  The  initial  intent  was  to  place  the  receiver  and  support  the  target  on  the  same 
platform.  However,  the  removal  of  the  target  for  the  background  noise 
measurement  caused  flexing  in  the  platform  and  produced  source-to-receiver 
movement  giving  incorrect  results  in  the  subtraction  routine.  Figure  9  illustrates 
the  effect  of  source-receiver  movement  on  "target-in"  and  "target-out"  waveforms. 
The  solid  trace  represents  the  target-in  waveform  and  the  dashed  line  represented 
the  target-out  waveform.  The  reflection  began  at  35  usee.  The  dashed  and  dotted 
line  is  the  the  result  from  the  subtraction  of  the  two  waveforms  and  it  was  clear  that 
in  this  case,  it  did  not  give  the  correct  results.  The  region  between  0  and  35  usee 
should  have  cancelled  to  zero.  The  non-zero  result  indicates  a  phase  change  caused 
by  relative  source-receiver  movement.  The  phase  change  subsequently  affected  the 
amplitude  of  the  reflected  signal  from  35  to  100  usee.  It  was  easy  to  see  that  any 
phase  change  brought  about  by  small  source-receiver  movement  would  cause 
erroneous  results.  To  eliminate  this  problem,  a  separate  platform  was  subsequently 
built  and  used  for  the  receiver. 

The  output  of  the  power  amplifier  was  also  checked  in  preliminary  work 
for  fluctuations.  An  arbitrary  source-receiver  geometry  was  selected  without  any 
target  in  the  tank.  Several  waveforms  of  the  direct  pulse  were  stored  at  various 
times  after  initial  start  up.  Amplitude  fluctuations  existed  shorty  after  start-up. 
Ideally,  successive  waveforms  subtracted  from  each  other  would  cancel  out  to  give  a 
zero  response.  If  amplitude  fluctuations  existed  a  false  signal  would  occur  and  this 
in  fact  did  occur  after  five  minutes  of  warm-up  as  Figure  10  illustrates.   Another 
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test  taken  at  15  minutes  after  start  up  produced  better  but  still  unsatisfactory  results. 
After  one  hour  of  warm-up  the  amplitude  fluctuations  subsided  and  a  zero  respose 
was  obtained. 
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IV.   FAR  FTELD  DETERMINATION 

The  measurement  of  the  far  field  of  several  objects  was  conducted.  An  equation 
that  gives  the  approximate  beginning  of  the  far  field  is  given  by  [Ref.  2:p.  188]: 

2 

rf=—  (3) 

AX 

L  =  the  longest  dimension  of  object 
X  =  wavelength 

It  is  always  best  to  determine  the  far  field  by  measuring  pressure  amplitude  as  a 
function  of  distance  from  the  object.  The  equation  for  a  diverging  spherical  wave, 
as  described  by  equation  2,  can  be  manipulated  by  taking  the  log  of  both  sides  to 
transform  it  into  a  linear  equation  with  a  slope  of  minus  one.  The  far  field  can  be 
determined  experimentally  when  the  observed  slope  of  the  curve  starts  to  deviate 
from  minus  one. 

The  far  field  of  the  F41  was  determined  to  begin  at  approximately  100  cm  and 
70  cm  for  176  kHz  and  70  kHz,  respectively.  The  far  field  for  the  sound  reflected 
from  the  styrofoam  cylinder  positioned  broadside  to  the  source  was  determined  to 
begin  at  40  cm  for  70  kHz.  The  three-inch  diameter  styrofoam  sphere's  far  field 
began  at  less  than  20  cm  at  70  kHz.  Figures  1 1  through  14  are  graphs  of  the  data 
used  to  determine  their  far  fields.  Lines  of  slope  equal  to  minus  one  are  drawn  on 
the  graphs. 
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V.   SPHERES 

A.  INTRODUCTION 

To  achieve  an  understanding  of  the  complexities  involved  in  analyzing 
scattering  effects  of  various  objects,  it  was  instructive  to  begin  with  a  very  simple 
and  symmetrical  object,  a  sphere.  Several  solid  metal  spheres  were  used.  Spheres 
were  aluminum  and  brass,  but,  the  exact  composition  was  not  known  or  investigated 
since  they  were  to  serve  only  as  preliminary  tools  in  studying  the  general  properties 
of  waveform  scattering. 

Since  metal  spheres  are  not  perfectly  rigid,  they  are  capable  of  absorbing 
energy  and  subsequently  reradiating  it.  Also,  they  possessed  resonances.  Many 
experiments  in  the  past  30  years  have  been  conducted  with  spheres  and  long 
cyHnders,  and  the  Resonance  Scattering  Theorem  (RST)  has  evolved  from  this  work 
[Ref.  5,6,7,8,11,12,13,14].  This  theorem  predicts  scattering  patterns  which  entail 
both  the  rigid-like  properties  and  the  resonance  properties  of  an  object.  It  was  not 
the  intent  to  examine  this  theorem  in  detail  but  only  to  observe  some  of  its 
predictions  experimentally. 

B.  MONOSTATIC  FORM  FUNCTION 

The  goal  of  replicating  some  of  the  data  in  the  literature  was  sought  to  establish 
confidence  in  the  experimentation  process.  Neubauer,  Vogt,  and  Dragonette  [Ref. 
7]  conducted  an  analysis  on  elastic  spheres  using  steady-state  signals  and  observed 
the  backscattered  signal  form  function  from  an  aluminum  sphere  as  a  function  of  ka. 
The  reflection  form  function  is  defined  as: 
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f=^ 


Ps 
Po 


(4) 


Ps  =  scattered  pressure  amplitude  measured  at  r 
Po  =  incident  pressure  amplitude  at  sphere  center 
a  =  sphere  radius 

Their  experimental  data  were  plotted  and  compared  against  theoretical  results 
and  good  agreement  occurred. 

Part  of  Neubauer's  experiment  was  repeated.  A  four-inch  diameter  aluminum 
sphere  was  ensonified  by  50  cycle  pulses  from  100  to  130  kHz,  and  the  backscattered 
pressure  amplitude  was  measured.  Time  discrimination  of  the  reflection  was  used. 
The  amplitude  of  the  reflected  pulse  was  taken  at  what  was  judged  to  be  the  steady- 
state  part  of  the  signal.  The  receiver  was  1.60  m  from  the  center  of  the  sphere, 
placing  the  receiver  well  within  the  far  field  of  the  sphere. 

The  wavelength  used  in  the  calcuation  of  ka  was  dependent  on  the  speed  of 
sound  in  the  water,  which  in  turn  was  affected  by  temperature,  pressure,  and 
salinity.  For  the  fresh  water  tanks,  salinity  was  ignored.  The  speed  of  sound  was 
based  on  an  empirical  formula  that  applied  to  distilled  water  [Ref.  l:p.  107]. 
Although  the  water  was  normal  tap  water,  the  distilled  water  formula  was  a  close 
enough  approximation  [Ref.  15].  Water  temperature  of  the  tank  was  a  consistent  17 
degrees  Celsius  throughout  the  entire  experimentation  phase.  At  this  temperature 
and  normal  atmospheric  pressure,  the  speed  of  sound  was  calculated  to  be  1473  m/s. 
Experimentally,  two  LC-IO's  were  placed  on  a  track  180  cm  apart  and  the 
successive  arrival  times  of  the  pulse  were  measured  via  the  Tek  2430  with  the 
cursor  controls.  The  speed  of  sound  thus  was  measured  to  be  1480  ±  20  m/s,  placing 
1473  m/s  well  within  the  experimental  uncertainty,  and  this  value  was  used  as  the 
speed  of  sound  in  water. 
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The  shadowing  effect  of  the  LC-10  in  front  of  the  F41  source  was  taken  into 
account  by  placing  a  dummy  LC-10  at  that  position  and  recording  incident  pressure 
at  the  sphere  center  location.  Additionally,  since  the  LC-10  was  not  perfectly 
omnidirectional,  the  support  rod  was  marked  to  maintain  a  constant  receiver  aspect 
with  the  target.  Placing  the  LC-10  in  front  of  the  F41  also  caused  additional 
reflections  to  occur.  The  direct  pulse  was  reflected  by  the  LC-10  back  to  the 
source's  face  and  again  reflected  there  back  to  the  LC-10.  More  "noise"  was  present 
when  measuring  the  monostatic  case  but  appropriate  location  of  the  receiver 
eliminated  any  interference  from  those  reflections. 

The  data  matched  Neubauer's  results  well,  and  the  form  function  as  a  function 
of  ka  is  plotted  in  Figure  15.  Data  from  Neubauer's  experiment  is  also  plotted  in 
Figure  15.  The  only  difference  between  results  was  a  small  shift  in  ka  of 
approximately  0.2-0.3,  not  surprising  since  the  compostion  of  the  sphere  was 
unknown.  The  minima  corresponds  to  resonant  frequencies  of  the  sphere,  where 
the  specularly  reflected  wave  is  interfering  with  the  resonantly  radiated  wave  from 
within  the  body. 

C.    WAVEFORM  EFFECTS 

The  effect  of  a  sphere's  resonance  was  examined  in  more  detail.  A  four-inch 
diamter  aluminum  sphere  was  ensonified  by  a  30  cycle  pulse  at,  above,  and  below  a 
resonant  frequency  of  82.7  kHz.  At  resonance,  the  ringing  of  the  sphere  is  clearly 
seen  in  Figure  16  when  the  direct  pulse  ends.  Pulse  length  was  approximately  367 
|is  long.  The  sphere's  stored  energy  is  simply  being  reradiated.  Figures  17  and  18 
show  the  waveform  response  2  kHz  above  and  below  resonance.  They  are 
decisively  different  and  clearly  illustrate  the  resonance  effect. 
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The  effect  of  sphere  composition  on  reflected  waveforms  was  observed.  A 
four-inch  diameter  brass  sphere  was  ensonified  by  the  same  frequencies  as  the 
aluminum  sphere.  Brass  has  a  characteristic  impedance  approximately  2.5  times 
that  of  aluminum.  Figures  19,  20,  and  21  show  the  reflected  waveforms  having 
frequency  dependence  that  are  markedly  different  from  the  aluminum  sphere. 

The  effect  of  pulse  length  was  observed  using  a  4-inch  aluminum  sphere  at  and 
off  a  resonant  frequency.  Figure  22,  at  resonance,  and  Figure  23,  off  resonance, 
show  reflected  wavforms  using  a  short  pulse  of  2  cycles.  They  appear  very  similar 
in  shape.  Figures  24  and  25  show  the  results  for  a  pulse  length  of  5  cycles.  Again, 
the  results  of  the  two  are  similar.  When  pulse  length  reaches  15  cycles,  the  effects 
of  the  resonance  begins  to  be  clearly  observed,  as  Figure  26  and  27  illustrate. 

For  comparison,  the  same  procedure  above  was  performed  on  a  four-inch 
diameter  brass  sphere.  Results  are  illustrated  in  Figures  28  through  33.  The  brass 
sphere  did  not  show  any  resonance  effect  as  determined  in  the  preceding 
composition  experiment  for  these  frequencies,  so  its  results  were  very  similar  to  the 
aluminum  sphere's  results  off  resonance.  Waveforms  were  practically  identical  for 
the  shorter  pulse  lengths  regardless  of  frequency  and  composition  implying  the 
reflections  are  purely  specular.  The  results  clearly  demonstrated  the  dependence  of 
pulse  length  on  reflected  waveforms. 

D.   SELECTION  OF  STYROFOAM 

Since  defining  target  strength  becomes  difficult  when  using  elastic  spheres, 
which  was  not  the  intent  of  this  work,  the  need  existed  to  do  one  of  two  things:  (1) 
work  at  a  frequency  that  gives  good  replication  of  the  input  pulse,  or  (2)  find  a 
target  that  gave  a  very  high  impedance  mismatch. 
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Figure  16.  Aluminum  Sphere  Reflection  at  Resonance,  82.7  kHz 
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Figure  17.  Aluminum  Sphere  Reflection,  84.7  kHz 


35 


•^0 


2L.1 


(.0 


-  .  2t:t     - 


-  .  4D 


..,-,  .,  ,.  ,  J   ,.  ,   ,    ,1,    , 

— — T • T  -r-l 1— 

1 — , — ^ — — 1 — 1 — 1 — , — r— . — f — , — 1 — ^ 

•                  1 

-  i 

■         i!        !     Ml 

-  l! 
1              ii 

1         1 

II     ' 

i    1      J    hi 

■ 

1 

1     ji  [i 

iWiMW^^^^^^^ 

:        Ijl             1    ' 

1             11    Lii'i  P* 
i'     1              if' 

■        ill     ■  t'  '  1 

'                1 

' 

■            i'      f         t 

II      Ii         1 

J 

:                 'h 

1 

■ 

'  '     1  '  1 

"    ■    ' 1          i     .     1     .     1     .     J     .     1     .     i     .     1    . 

1  .  1   ,  )   .  t  ,   k  .   1  ,   1  , 

0 


!0£»  4  ©fci  St?Pt 

(micro sscands  i 

Pel  ail  ve     T  i  rne 


I  00 


Figure  18.  Aluminum  Sphere  Reflection,  80.7  kHz 
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Figure  19.  Brass  Sphere  Reflection,  80.7  kHz 
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Figure  20.  Brass  Sphere  Reflection,  82.7  kHz 


04 


0S 


QQ     U,. 


-  .  k5^ 


-.  E?4 


-I — ' — I — ' — r 


!  li 


ill 


m 


11' 


"ij     ilM 


I       l|  ill 

iii 


'>ii 


'ii 

i:  i.'l  I 

ijil    tl 

■  I! 


M 


if 


t),.i\fi)f,^'f,?,!llfi|fif,i„,,; 
|nvyiiiii.'v'A'!!Vii;^^H'vV^'^ 


^r'i'ipi' 


_i . I , ; , '     ■     '     ■     I     ■     I     ■     I 


J ■      I      ■ I 1 L 


2©G  4GG 

f  (Ti  1  c  r  o  s  e  c  a  n  6  s  I 

Re  I  a"t.  1  ve     T  itrie 


60© 


Figure  21.  Brass  Sphere  Reflection,  84.7  kHz 
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Figure  22.  Aluminum  Sphere  Reflection,  82.7  kHz 
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Figure  23.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  24.  Aluminum  Sphere  Reflection,  82.7  kHz 
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Figure  25.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  26.  Aluminum  Sphere  Reflection,  82.7  kHz 
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Figure  27.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  28.  Brass  Sphere  Reflection,  82.7  kHz 
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Figure  29.  Brass  Sphere  Renection,  84.7  kHz 
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Figure  30.  Brass  Sphere  Reflection,  82.7  kHz 
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Figure  31.  Brass  Sphere  Renection,  84.7  kHz 
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Figure  32.  Brass  Sphere  Reflection,  82.7  kHz 
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Figure  33.  Brass  Sphere  Renection,  84.7  kHz 
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The  most  important  factor  in  achieving  a  good  reflection  is  the  impedance 
mismatch  between  the  fluid  and  object.  Experimental  work  done  with  metal  spheres 
by  Hamptom  [Ref.  5]  shows  the  effect  of  compostion  using  several  metals.  A 
pressure  wave  in  water  will  see  little  resistance  to  a  hard  rubber  sphere  because  its 
characteristic  impedance  is  near  that  of  the  water.  An  object  that  gives  a  high 
impedance  mismatch  is  an  air  bubble.  The  impedance  ratio  of  water  to  air  is 
roughly  3000.  For  the  case  of  an  air  bubble,  the  boundary  condition  between  the 
water  and  air  is  termed  "pressure  release",  since  the  pressure  transmitted  into  the  air 
is  nearly  zero  [Ref.  2:p.  127].  The  difference  between  a  pressure-release  sphere  and 
a  rigid  sphere  is  best  seen  graphically  in  Figure  34  [adapted  from  Ref.  15:  p.  359]. 
In  the  ka  region  used  in  the  experiments,  both  the  rigid  and  pressure-release  form 
functions  will  be  nearly  identical  and  equal  to  one. 
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Figure  34.  Backscattered  Form  Function  of  a  Rigid  and  Pressure 

Release  Sphere 

A  styrofoam  ball  was  used  to  approximate  an  air  bubble  not  near  resonance. 
The  reflected   waveforms   were  generally   identical  to   the   input  pulse   at   all 
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frequencies.  This  made  styrofoam  an  ideal  target  for  this  work.  Figures  35 
through  40  show  reflections  from  a  three-inch  diameter  styrofoam  sphere  using 
various  pulse  lengths  at  100  and  106  kHz. 

E.  BISTATIC  SCATTERING 

A  360  degree  pattern  of  the  pressure  amplitude  scattered  from  a  styrofoam 
sphere  was  obtained  to  compare  the  results  with  theoretical  predictions  for  a 
perfectly  reflecting  sphere.  Since  a  sphere  is  perfectly  symmetrical,  it  was  not 
necessary  to  take  the  entire  360  degrees  of  data  but  only  1 80  degrees  of  it. 

Theoretical  equations  have  been  derived  and  published  by  several  authors  that 
predict  scattering  of  acoustic  plane  waves  from  a  sphere  [Ref.  16,17].  They  tend  to 
be  quite  complicated  functions  involving  sums  of  Hankel  functions.  Figure  41  was 
adapted  from  Medwin  and  Clay  [Ref.  16:pp.  192-193]  and  shows  computed 
scattering  patterns  in  relative  amplitude  for  a  perfectly  reflecting  sphere  in  water 
for  several  ka  values.  The  incident  plane  wave  is  arriving  from  the  right.  The 
patterns  show  that  for  increasing  ka  the  major  forward  lobe  gets  larger  in  amplitude 
and  narrower  in  width.  More  structure,  i.e.,  secondary  lobes,  develop  as  ka 
increases. 

Bistatic  data  were  obtained  for  a  ka  =10  using  a  three-inch  diameter  styrofoam 
sphere  and  are  shown  in  Figure  42.  The  incident  wave  arrives  from  the  bottom. 
Frequency  was  61.5  kHz  and  an  eight  cycle  pulse  was  used.  The  graphed  data 
revealed  some  structure  in  the  forward  scattering  region  and  a  large  forward 
scattering  lobe.  Comparison  with  the  theoretical  prediction  for  a  rigid  sphere 
(dashed  line  in  Figure  42)  show  that  the  salient  features  of  the  scattered  pattern 
were  duplicated. 
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Figure  35.  Styrofoam  Sphere  Reflection,  100  kHz 
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Figure  36.  Styrofoam  Sphere  Reflection,  106  kHz 
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Figure  37.  Styrofoam  Sphere  Renection,  100  kHz 
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Figure  38.  Styrofoam  Sphere  Reflection,  106  kHz 
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Figure  39.  Stvrofoam  Sphere  Reflection,  100  kHz 
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Figure  40.  Styrofoam  Sphere  Reflection,  106  kHz 
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Figure  41.  Theoretical  Scattering  Patterns  for 

Rigid  Spheres  in  Water 
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Figure  42.  Bistatic  Form  Function  of  a  Stvrofoam  Sphere 

ka  =  10 
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VI.  BTSTATIC  SCATTERTNG  OF  A  STYROFOAM  CYLINDER 

A  styrofoam  cylinder  of  diameter  2.5  cm  and  length  20  cm  with  flat  ends  was 
used  for  bistatic  scattering  measurements.  The  particular  diameter  was  chosen 
because  the  styrofoam  was  commercially  available  in  that  size.  The  length  was 
selected  using  two  criteria:  (1)  the  length  was  small  enough  so  that  the  receiver 
could  be  placed  in  the  far  field  of  the  sound  reflected  from  the  cylinder  and  (2)  it 
was  large  enough  to  clearly  present  a  cylindrical  shape. 

Since  the  F41  was  not  omnidirectional,  there  was  concern  of  approximating  a 
plane  wave  at  the  target.  The  lower  the  frequency,  the  better  the  pressure 
distribution  across  the  target  because  of  the  larger  major  lobe,  but  the  wavelength- 
to-target  length  ratio  would  be  less  and  the  signal-to-noise  ratio  also  became  poorer 
at  lower  frequencies.  A  frequency  of  70  kHz  was  chosen  at  a  source-to-target 
distance  of  125  cm.  Figure  43  shows  the  incident  pressure  distribution  as  a  function 
of  distance  off  the  axis  at  a  distance  of  125  cm.  The  lower  curve  is  for  the 
mono  static  case  where  the  receiver  shadowed  some  of  the  incident  energy  with  a 
receiver-to-target  center  distance  of  50  cm.  A  two  percent  pressure  amplitude 
difference  existed  for  the  shadowed  case.  An  eight  percent  decrease  in  pressure 
ampUtude  existed  when  not  shadowed. 

The  cylinder  was  mounted  to  a  1/8-inch  diameter  metal  tube  which  had  a  2  X 
1/2  X  1/16-inch  flat  plate  epoxied  to  it  to  provide  a  means  of  rotating  the  cylinder. 
The  flat  plate  was  forced  into  the  interior  of  the  cylinder  and  provided  a  snug  fit. 
The  effect  of  the  plate  on  the  reflections  was  assumed  to  be  neglible  since  it  was 
completely  surrounded  by  the  styrofoam.  The  tube  was  as  small  as  possible  to 
minimize  reflections.  A  string  with  weights  was  attached  to  keep  the  st>Tofoam 
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Figure  43.  Pressure  Distribution  at  70  Khz 

cylinder  submerged.  The  thin  tube  was  soldered  into  a  3/8-nch  diameter  tube 
which  supported  the  entire  arrangement.  The  cylinder  was  level  before  entering  the 
water  and  appeared  to  remain  that  way  when  in  the  water.  Figure  44  illustrates  the 
support  setup . 

Data  were  recorded  for  nine  bistatic  angles.  Time  discrimination  of  the 
reflected  pulse  was  used  whenever  possible.  The  source-to-receiver  distance  was 
125  +  2  cm  in  all  cases  with  the  F41  being  driven  at  177  V  rms.  The  Ithaco 
preamplifier  setting  was  20X  with  a  3-300  kHz  bandwidth.  Receiver-to-cylinder 
center  distance  varied,  but  the  receiver  was  always  in  the  far  field  of  the  cylinder. 
All  data  were  extrapolated  to  one  meter  from  the  target  center  for  comparison.  A 
20  cycle  pulse  was  used. 

The  cylinder  was  initially  aligned  broadside  to  the  source  by  eye  to  within  an 
uncertainty  of  about  +  4  degrees.  Assuming  symmetry,  data  were  taken  for  180 
degrees  of  rotation  and  several  points  were  taken  beyond  t  hat  limit.    Figure  45 
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Figure  44.     Support  Platform  Setup 
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illustrates  the  angles  measured  in  the  experiments.  The  varied  angle,  A,  is  the  angle 
made  between  the  cylinder's  normal  and  the  line  connecting  the  center  of  the  target 
to  the  source  axis.  Angle  B  is  the  bistatic  angle. 


Target 


Receiver 


Figure  45.  Measured  Angles  in  Experiment 

The  results  of  the  experiment  are  shown  in  Figures  46  to  55.  In  most  cases, 
strong  peaks  were  observed  when  angle  A  was  one-half  the  bistatic  angle  and 
secondary  peaks  with  smaller  amplitude  occurred  90  degrees  from  that.  This 
geometry  seems  to  indicate  specular  reflection  off  the  long  cylinder's  surface  and  its 
flat  ends  respectively.  As  the  bistatic  angle  increases  the  secondary  peak  broadens 
and  eventually  disappears.  Well  into  the  forward  scattering  region,  greater  than 
150  degrees  of  bistatic  angle,  the  data  became  less  predictable  The  distinction  of  the 
primary  and  secondary  peaks  became  less  clear  with  additional  peaks  showing  up. 
At  a  bistatic  angle  of  1 80  degrees  no  distinctive  peaks  were  present  and  amplitude 
became  very  large  in  comparison  to  other  bistatic  angles. 
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VII.   CONCLUSION  AND  RECOMMENDATIONS 

A  data  acquisition  system  using  the  Tektronix  2430  digital  oscilloscope  was 
implemented  to  extract  scattered  waveforms  in  a  noise-filled  water  environment. 
The  system  was  verified  to  work  and  reflected  waveforms  were  collected. 
Preliminary  work  revealed  that  the  relative  movement  of  the  source  and  receiver 
between  collection  of  waveforms  was  critical  in  obtaining  accurate  data  using  a 
digital  sampling  system  with  a  subtraction  routine.  Data  collection  was  slow  and 
tedious. 

Several  aspects  influencing  scattered  waveforms,  including  target  composition, 
frequency,  and  pulse  length  were  observed  and  recorded.  The  effects  of  these  on 
the  waveforms  of  the  reflected  pulse  produce  difficutlies  in  measuring  the  intensity 
of  pulses  with  waveforms  not  identical  of  the  incident  pulse. 

Styrofoam  targets  were  used  to  circumvent  the  problem  of  waveform 
dependence  on  frequency,  pulse  length,  and  composition.  Bistatic  measurements 
were  taken  using  a  sphere  and  a  cylinder  as  targets.  The  data  from  the  sphere 
compared  well  with  theoretical  predictions  for  a  perfectly  reflecting  sphere. 
Readily  available  data  were  not  found  for  comparison  with  the  results  obtained  for 
the  cylinder.  The  goal  of  matching  wavelength-to-target  length  ratios  observed  in 
oceanic  conditions  was  not  done  in  this  work.  The  USRD  type  E27  transducer 
would  be  an  excellent  source  for  this  application  since  it  has  a  frequency  range  100 
to  500  kHz.  Its  small  size  makes  its  beam  widths  very  large  which  is  beneficial  when 
receiver-target  distance  is  limited  by  tank  dimensions.  The  E27  is  also  very 
lightweight  and  easily  moved.  The  effect  of  the  jitter  would  have  to  be  evaluated  at 
the  higher  frequencies. 
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The  bistatic  radar  cross-section  theorem  could  not  be  compared  to  the  data 
collected  in  this  report.  The  wavelengths  used  in  the  experiments  are  not  close  to 
the  "vanishing  limit"  as  stated  in  the  theorem.  One  common  trend  between  the 
bistatic  data  and  the  theorem  was  that  the  data  became  less  predictable  as  the  bistatic 
angle  approached  180  degrees. 

For  future  work,  a  physical  structure  must  be  made  to  facilitate  positioning  of 
the  source,  target  and  receiver  in  three  dimensions  so  that  data  collection  can  be 
completely  automated  and  accuracy  improved.  The  rotation  of  supports  should 
have  a  digital  device  to  allow  control  by  the  computer.  The  design  and  construction 
of  such  a  system  would  not  be  a  trivial  project  and  would  require  a  great  amount  of 
ingenuity  and  experience  with  electronic  devices.  The  payoff  of  having  such  a 
system  would  be  great  to  students  studying  underwater  acoustics  and  in  particular  to 
anyone  studying  scattering  effects.  Once  a  positioning  system  is  built  and 
completely  automated,  data  collection  will  become  simpler  and  more  time  can  be 
devoted  to  analyzing  it. 
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APPENDIX  A 

COMPUTER  LISTING  OF  "TEX" 


4<d 

50 

60 

70 

80 

90 

100 

1  10 

1Z0 

130 

U0 

150 

160 

170 

180 

190 

Z00 

Z10 

2Z0 

Z40 

Z50 

Z60 

Z70 

Z80 

Z90 

300 

310 

3Z0 

330 

340 

350 

360 

370 

380 

330 

400 

410 

4Z0 

430 

440 

450 

460 

470 

480 

490 


ti  II  ii  tt  li  II  ii  II  11  II  '*  "  "  *'  **  **  **  * 

nniTirnTririin  uiiitPTritiiTTT 

II  11  »i  rr  Ft  TT  II  II  tt^t^t^^^^t  II  jf  ^r^^'T 


PROGRAM   TEX 


*      M  t      *fc      **       M   ±      M   t      *■       **      **      *fc       fcJ       *|       **      *■      **      *1       aj      **      **       ■*      ji       |J       t    L      ■*      **      *J       **      *J       *A      JA      ±  ±      *|       ^1      i|      ■*       **       ■* 


DRTft  AQUISITION  FOR  BI-STATIC  TARGET  STRENGTHS 
USING  THE  TEKTRONIX  Z430  DIGITAL  OSCILLOSCOPE 

THIS  PROGRAM  TRANSFERS  WAVEFORM  DATA  FROM  THE  TEX  Z430 
TO  THE  COMPUTER  UHERE  IT  IS  CONVERTED  TO  ACTUAL  VOLTAGE 
VALUES  AND  STORED  ON  A  FLEXIBLE  DISK.   INPUT  FROM  THE 
USER  REQUIRES  FILE  NAME,  DESCRIPTION.  AND  CONTINUATION 
RESPONSES.   PREVIOUS  SETUP  OF  THE  TEX  Z430  IS  REQUIRED. 
A  THOROUGH  KNOWLEDGE  OF  THE  GPIB  CODES  IS  NECESSARY.  SEE 
THE  OPERATOR'S  MANUAL  AND  THE  INTERFACING  MANUAL  FOR 
THE  TEX  2430. 

■  *  ■  *   ■  *  *  ■  *  *   fc  ■  m  t      it      *  t  \  m     i  A     ■*  a  a      *■  a.  k     ±  A     j*  *m   ■■  **  at     a  m      at  J  fc  ■  *   *  ■  A  %.    A  a      at     t  a     m  t      t  a      t  t      a  a     JA  t  a      t  t      at      *  ■   ■  ■   *  a     t  i  j*  ■■   ■*  A  a     j«  A  a      at      a  a     t  a      at      t  a      a.  a  a      ■  ■   *  .   ■  *   A  J  ^^  ^^..^^^^L^^ 

HnTtitTrTrTTn  TtiTiiniinTi  u  itannnira  Itlt^ 

^T  Tl    fW  ^TTT^TTT^^^Tt^T^tTT^XTT^T^t^t^T^T^^^^^T^T^TtT^T^TT^tt^X^XTT  TF  TT  ?I  ti    ti  TT^T^T^tTX  T?  ti  ^T^t  ti  ^  I  TT  FI    TT   FT    TI    II   TTFl    II  TT  II    If   II    II  TTTT 

INITIALIZATION 

*  *  at       ■■   J*   A*   ■.   ■■   AA  at      At      t  t      a    a       at      h.   aA  a   a       t   t       >*   ■■   ■■   ■*   a1  at       it   Al  t    a      A*   IM   it   Al   AA   *fc   ■*   ■*   ■*   *■  »    k      a    a       ■*   **   AA  t  t       at       AAA*   Al   ■■   A*   A*   lA   t*   Jt   tA   *A   ■*   AA   lA   ■*   AA   AA   MA   ■■   AK   ■* 

HTnTauiiiriTTrii  ii  n  nn  iiim  n  iiniiirairiririiiTnTrTiTiTTairnitiTiriTnTia      ii  n  uiiiTTrn  iriiiTiriiiriTagirTiir 

^T  tT^t^T^T  tTTTTT  II    fTTTTT    II    IITITIIT  TVtT  T^  TT   IT  T  T  T  T  TT  TTTtTT  FT   IT  TT^TTT  FT  TF  TT  ^T  TT  TT^^T^W^T^^  FI    II    II    II  ^tTTTT  T  I   TT  ^T  TT  ^T  Tt^t  TT  T^^T  1 1    II  TT  TT 

PRINTER  IS  CRT 

OPTION  BASE  0 

DIM  Vdata$[ 15] ,U$I40] !NAME  OF  DATA  FILE  AND  DESCRIPTION 

DIM  A$[Z00] ,B$[5000] ,D$[5] IPRE-UAVEFORM  AND  WAVEFORM  STRINGS 

DIM  V( 1023), Tin( 10Z3) ,Volt( 10Z3) 'DATA  ARRAYS 

DIM  X$[ 10]  ,Y5$[  10]  ,Vs$[  10] ,Ts$[  10] 

DIM  An$[  1  ] 

AA   .1   lA   A.   .1   1.   11   .1   .1   lA   AA   ..   ».   1|   AJ   AA   lA   AA   tl   .1   ..   Al   1.   A.   H   JA   11  tt      t.       »■   ..   Al   lA   11   Al   AA   AA   ll   .1   .■   ll   1.   I.   ..   Al   II   J.   .1   II   ll   11   II   ll   ll   11   II   11   ll   A  I   11   ■  ■   1  1   *  ■   '  ■ 

iiiiiiiriin  unTTirn  nil  TT  n  iiiiiriTnniinTriTiTTTiT  niT»  uitTraTrnTi  itniinirirTtTTn  n  iiiriririrpii  irirnTrTTiiTTTrn 

TT^X  TT  ^T^T  II    II    II    II    II    II  TT  TT  TT  ^T  ^T^T^T^  TT  TT  ^T  II    II  TT  TT  TT  TT  TT  TT  ^T  TT^T  ^T  TT^T  TT  TT  TT  IT  TF  FI   II    II    TT  TT  TT  ^TTT^T^  II    II    II    II    II    II  TT  TI   TT  ^T  I  F  IT    F  I    TI    II 

TEXTRONIX    SETUP 

AA   Al   Al   11   Al   *1   AA   lA   A*   Al   Al   11   11   lA   ll   ll   lA  t  a       at      At       11   lAAl   ll  at      AA   lA   11   11   ll   AA   ■..!   AB   AA   AA   11   lA   ■!  M   t      AA   Al   Al   Al   AA   Al   AA   11   Al   AA   AA   11   Al   Aj   lA   lA   AA   lA   Al   Al   AA   *1   lA   AA  MM 

HTTTTTTTT  tt  H  IT  rrTTTrTTTTpnTI  ITTrTf  n  TTTTtTTTTinnTITI  TTTTTTTTTI  TTTTtTtT  TT  TTTf  TT  TT  TTtT  TT  TTiT  IT  a  IT  II  II  II  TTITmrTTim  TT  TI  TT 

TTTT  TI    TI    FI    FFTF  TTTT  FI   TI  ^T^^  FI    II    FI   II    II    FI    II  ^tTT^T  TTTTTTTI    IIIF  TTXT^T^TTT^W    II    II  ^ttT  II    II    II  ^T  II   TT    FTTITI   TI    TI    II   II    FI    II   FI    II    IF    FI    TF   FI    ITTF  TTTT 

OUTPUT  704; "PATH  OFF" ! SHORTENS  THE  DATA  STRING  SENT  BY  THE  Z430 

Al       lA       AA       11       lA      .   a       Al       «■       Al       11       11      Al       lA      Al      A.       .1       AA       Al       lA      AA      AA       AA      Al       AA       AA       Al       lA      AA      AA      AA       Al       Al       Al       11       AA       lA       11       11       AA       Al       AA       lA       lA      AA      lA      AA       AA      11       11      Al      11      11      AA     Jl      Al       11      AA      AA      AA      lA      U      Al      Al     JA 

itiTiiiiititTiTr  itnu  TiiTirTiii  itttiiiriinTtirii  11  iiTtiiirii  uTrTtij  irtrit  utiiiii  HitiTiiTi  iTii  iiirirtriririiiiiTiriiirTrTiTT 

TT  TT  "^F^W^W  ^T  ^T  TT  TT^T"^T  ^T  TT  ^T  ^T^T  IT  Tf  TT^^^I    FI    IF  TT  ^TTT^^^T^T^T  TT  ^T  II    FT   FF   FF    FF    II  ^T  IIFFTF    II    IIITTI  ^TTT^T    II   TT  TT  II    II    II    II    l|    II     II    II    IF  TT  TI   TI 

INPUT  FROM  USER 

■  1   lA   11   11   11   11  at       11   AA   lA   Al   11   11   11   11   11   11   11   Al   11   ll   AA   Al   Al   AA   ll   11   AA   AA   ll   AA   11   11   11   11   11   11   11   Al   11   Al   AA   Al   11   11   11   Al   ■*   ■■   lAAl   Al   1-1   AA   lA   lA   AA   Al   11   lA   Al   Al   Al   1.A 

iriiiinHiiiriiTiitii  TrTrTTHiiiT  u  itiihb 

TT  II    FI  ^tTT^T^T^T^T^t^TTTTF^T^T^T^T^^^^^^^X^T^T^^T^T^T^T^t^TTT^T^^X^T^X^T^W^T^T^X^T^T^T  TTTF   TI    II   II    II    ITTTVT  TT^^TT^TTTTT  TI    FT   TT  TT    FT    FF    FT 

INPUT  "NAME  OF  THE  DATA  FILE?" . VdataS 

INPUT  "DOES  FILE  ALREADY  EXIST?( Y=YES ) " , An$ 

IF  An$="Y"  THEN  PURGE  VdataS 

INPUT  "INPUT  FILE  DESCRIPTI0N7( 40  CHARACTERS  MAX)".U$ 

BEEP 

PRINT  "ALLOWING  DATA  DISK  INSERTION.  PRESS  CONTINUE  TO  RESUME." 

PAUSE 
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500 
510 
5Z0 

530 
540 
550 
5G0 
570 
580 
590 
600 
B10 
GZ0 
B30 
640 
650 
660 
670 
680 
690 
700 
710 
7Z0 
730 
740 
750 
760 
770 
780 
790 
800 
810 
8Z0 
830 
840 
850 
860 
870 
880 
890 
900 
910 
9Z0 
9  30 
940 
950 
960 
970 
980 
990 
1000 


CREATE    BOAT    Vdata$ . 1 , 1 7000 
ASSIGN   eOfile    TO   VdataS 
PRINT 


».       ««       »t      H      ->      ».       .«       .»       ..       ..       ..       J.       «.       ..       ^1       It       1.       ■■       »m       .»       J.      Jl      AM       »»       .t      »»       »»      «»      J.      »»      A«      »J       ■»       »»      J»       ■■       ■.       -»       ».       »^      «1       »J       U      J»       1»       ■»       ..       ..       ..       ..       »»       ..      »«      J.       1»       «1       »j       .1       tt       ..       .»       ..       »«       ». 

TiiTii  nitiiirtriitrirtwttttiiiniiirirn 

OBTAIN   WAVEFORM 

*-*  **  **  **  **  ■■   ■*  **   ■*  **  fc*  **   fc*  *»  **  **  -*   **  ■*  **  *J   ■*  ■*   ■■  **  ■*  **  ■■  **  **  ■■   ■*  *i   a*  ■**!   ■*   ■*   **   *-  **   ■*   **  ■*  *■  fcM  jfc  ■*   at  J*  *■   t*   ■*  A*  *j   *■   *■   fct  *fc  *fc  t*  ■■  a*   ** 

iLii  iTit  irit  iriiiiirii  II IX  iriuuCitiuuuut  in 


PRINT    "OBTAIN   DESIRED   UAVEFORM   ONTO   TEXTRONIX" 

PRINT    "PRESS    CONTINUE    TO   RESUME" 

PAUSE 


■   ■       *■       **      J*     *1      **     Aa.     AA      AJ      a*     J*      A-*     J-l     J-A    -fc^    AA     A.*.    A  A.    J>A     A*     A-*.    .*-■■    -A*      J-*.    .AJ.   AA      A  1      AA    AA     AJ      *J.     ±  ^      *-^     A*.     AA     .**      **.     AA     A  A      AJ      A  J.    J.  J      AJ.     A  A     A-*.     A  A     AA     >  *.     A  A      fcJ     AJ      A*       BA     AJ      aA      AJ      A*      ^J     AA     AA      A*       fcA     JA     AA 

•trTrTr  TT  TT  Tl  TT  TT  Tl  H  t"  TTTTn  IIII  IT  II  IT  ITTmilTriT  11  II  ilirn  U  tTH  TI  TTTf  FT  TT  H  TTtT  TT  Tf  TI  TT  TI  TT  If  II  n  TTTTITtI  11  II  ll  O  TTTriT  Tl  ll  H 

^T^X^T^^^^r^^^T^^r^^^^r  f  f  If  ff  ¥▼  II  11  II  ff  TFff  FT  II  If  II  ^T^^r^T  fi  IP  II  11  II  II  ^^^^r^^^^r^^r^^r^rTI  11  ^X^^^^^^T^^r^X^^^X^X^^^^^^r^T^T^^r^^^^^T^^r^f  1 1  ^x 

RECEIVING   PRE-UAVEFORM    DATA   AND   THE    UAVEFORM 

■    *       **       aj       *J       tJ       MJ       *1       *1       A*       t*       *t       ^fc*i       **        tt       a*       *!**       *■       t*       *■       ■■       **       **       **       tJaJ       AAtlAAAAAA       AaAI       i^AAAJ       AA^aAAAAAAAAAJAAAAAA     A-A      AA     A.  A      A  A      A  A 

TTTrTTTrTT  Tl  IIII  uTIliTI  n  n  mrH  "  ^^  "  TTTTITTT  TTTTTITTTTTTTrTTTTTTTTTTTTTTTTTTTTTT  TTTTTT  n  ilTTTIII  ITlT 

OUTPUT   704; "UFM?    XIN"!TIME    INCREMENT   BETWEEN   POINTS 

ENTER   704; X$ 

OUTPUT   704;"UFM7    YMU" I  VOLTAGE    DIGITIZING    INCREMENT 

ENTER   704;Ys$ 

OUTPUT   704; "CHI  7    VOL" ! VOLTS/DI V 

ENTER  704; Vs$ 

OUTPUT   704;"HOR7    BSE"!B   TIME    SCALE 

ENTER  704;Ts$ 

! 

I 

OUTPUT    704;"CURV?" 

ENTER   704 ;B$ 

am     **      ■*     afc     at      *■      ■*      >■      at      ||     j|     ja      *t      jj      **     *j      *■      t*      ia     **      *■     at     aj      **      **     J*      *■      ■*      **     *■      ■*      ■»      fc*      *i      ■*     *■      *■      i*      **     *J     A*     Ji      |*      i|     A|      tl      11      *J     *a     J>      **     A*     JJ     J*     J*      ■*      ai      *■    ^^^^^^^^^^^^ 

trtrttittrgirtTttti  nn  rtnirniiTrirtitrTiuTrHtrtiltiriiiiliiLtt 

II    r  f  ^T^T  ^T^T^T  TT  II    II   TT  TT  TTTT^T^T^T  fl    II    II    rr  IT   IT   TT  ^TTT  TITFff   fl    II  ^TTT  fl    II    II    II   Tl  TT    II    If   II    II    IT   fl    II  ^TtT  IFTI   fl    |l    II   fl   II    II  TT  II    II    II   II   II    II  TT 

SEPARATION   OF    ASCII    STRING    INTO   DATA   POINTS 

ai      i*      *■      a*a*      ■*      **a*      at      aaaaaa       i*      A*      ■■      a*      aa      >*      *>      **■■      itai      *aA*aa      !■      !**■      k  ^      **      >i      *a      atiaai      i*      aaaaaaaa      ■*      aaa*      >*ai      >iA.t    ^a      !■      a*      aa     aa     ai     ai    ^^^^^^^^^^^^^^^^^..^^ 

itTinn  tTTiTTTtit B  u  TtTritJtPtrtTlinna  a  U  unnirtilTitTtlTTTTttrTru  a  iTn  It  Tl  TTH  tTTT  tt 

TT  TT  TT  TT  TT  TT  T I  T  I  TT  ^T^T    II    II   fl   II   II    II    fl    Tf  TT  ^T  ^T  If    II    II    II  ^TTT^T  ^T  TT  ^T  TT  fT  Tl   II    II  TT  ^T  TT  II    II    II  TT  TT  TT  IT  T  T   II    II  TT  II    II    II    II    II    11    fl    IT   II   fl  TT  TT  TT 

N=LEN(B$) 
D$="" 
lnc=0 
I 

FOR    1=1    TO   N 

IF    B$[I,I]=","    THEN 

V( Inc)=VAL(D$> 

Inc=Inc+l 

D$="" 

GOTO   950  ** 

END    IF 

D$=D$&B$[I.I] 

NEXT    I 


■  a  aa  aa  aa  ^  aa  aa  aa   >■   ja  ■■  %  t      aa   aa  mm     aa  aa  aa   aa  aa   aa  lA  a*   aa  aa   aa   aa   aa  Mt     aa  ta  ia  ja   la   11  aa   ■■   aa  ai   aa  aa  aa  aa   aa  aa  aa  aa  aA  aA  aa  AI  Ai  AA  Aa  aa  aa  Aa  Aa  Aa  Aa  aa  aA  AA  A.A 

TTTT  TTtTTI  Tl  TT  TT  TT  TrTTT!  TT  TT  TT  TTtTTTTTttTT  IT  Tt  Tl  Tl  fTTTTt  tt  TTII  tT  Tl  ttittTtT  TT  TrTrTT  TT  TT  Tl  TTIIIITrnTTII  ITTTII II  TITT TT llTI  II  \l  H  II 

II    II    II    II    II    II    II    II   fl   ff   II   Tl   Tf   If   II    II    II    II  ^T  1 1  TT   f  T   f  I   Tf   TT  T?  TT  TTTT  II    II   II    fl  ^T^T'^T  FT   TT  ^T  If   Tl    fl    fl    II    II   If   fTTf  TTTT  TTTT    IT   II    II    II   Tl   Tf  TF   Ff   ff    II    II    II 

CONVERSION   TO   REAL-TIME    DATA 

at  al  aa  aa  .A  aa  Aa  aa   aa  aa  aa  ai  aa   aa   a.  a.  ..   .1   .a   ■■   aa  aa  aa   aa   at   ti   aa   .■   aa  at  aa  a  a   at   aa  aa  aa  aa   at  aa  at  aa  ..   aa   aa  aa   ta  aa  It   a.   aa   aa  at   ta  tt  aa  aa   aa   aa   It  At  al   aa  AAAA 

nirirnu  mi  itTTTrinritiiiin  ii  ttti  Ttitit  ii  p  irit  nTTanniiiiitiririiiiTtiiTTiriiirn  uii  ixitiiii  nTririinirii  uTni  n  mr 

TT  TT  TT  I T  TT  IfTI    IF   fl    If   If  TTTT  f  I   IF   f  I    f  f  TT  II    fl   fl    II    II    II  TTTT  ff    II   11   ff  ^TTT^TTT  TT  Tl  TTTTTF  TT  II    II  TTITTTTT  II   II    II    II    II   fl    II    II   Fl    II    Fl  TT  II    II    II    If    II    II 
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1010 

10Z0 

1030 
1040 
1050 
1060 
1070 
1080 
1090 
1  100 
1  110 
1  120 
1  130 
1  140 
1  150 
1  160 
1  170 
1  180 
1  130 
1Z00 
1Z10 
12Z0 
1Z30 
1Z40 
1Z50 
1ZG0 
1Z70 
1Z80 


X$=TIME  INCREMENT  BETWEEN  POINTS   DIV/50 

YS$=VOLTS  BETWEEN  POINTS   DIV/Z5 

VS$=V0LT5/DIV 

TS$=B  TRIGGER  TIME  SCfiLE 


Del-tin=VAL(X$) 
Delvol=VAL(  Ys$) 
V5cale=VfiL(  Vs$) 
Tscale=VfiL(Ts$) 

I 

OUTPUT  @Df ile;U$.Tscale,Vscale 

FOR  1=0  TO  10Z3 

Timd  )  =  I»DeltiPi 

Volt( I )  =  V(  I >»Delvol 

NEXT   I 

! 

PRINT  "URITING  TO  FILE" 

OUTPUT  @Df ile;Tin(«) ,Volt(  *) 

RS5IGN  eOfile  TO  • 

BEEP 

I 

In    At   tk  ki   tl    II    '*    **    **    *'    **    *'    "    '*    **    *'    '*    '*    '*    ''    *'    '*    *'    '^    '*    ''    **    '*    "    '*    "    '*    '*    ''    '*    **    '*    *'    *'    *'    *'    '*    *'    *'    '*    *'    *'    *'    '*    '*    "    '*    "    *'    ^  '    '  •    **    "    **    '*    ''    *'    '*    " 

! 

INPUT  "DO  YOU  URNT  ANOTHER  UftVEFORM^  Y=YES".Bn$ 

IF  Rn$="Y"  THEN  GOTO  430 

I 

PRINT  "PROGRRM  TERMINATED" 
END 
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10 

Z0 
30 
40 
50 
60 
70 
80 
30 
100 

1  10 
120 
130 
140 
150 
160 
170 
180 
190 
Z00 
Z10 
ZZ0 
Z30 
240 
250 
260 
270 

2  80 
230 
300 
310 
320 
330 
340 
350 
360 
370 
380 
330 
400 
410 
420 
430 
440 
450 


APPENDIX  B 

COMPUTER  LISTING  OF  "PLOT" 

PROGRAM:    PLOT 

THIS   PROGRAM   TAKES   TUO   FORMATTED  BOAT   FILES   AND   SUBTRACTS 
ONE    FROM    THE    OTHER.       IT   UA5   DESIGNED   FOR   USE    USING   WAVEFORMS 
OBTAINED   ON    THE    TEKTRONIX   2430   DIGITAL    OSCOPE.       DATA    IS    TAKEN 
FROM   THE    FILES    IN    THE    FOLLOWING   ORDER. 

W$   =    FILE    DESCRIPTION   STRING 

TSCALE    =    TIME    SCALE    OF    THE    SCOPE    (SEC/DIV) 

VSCALE    =    VOLTAGE    SCALE    OF    THE    SCOPE    (VOLT/DIV) 

T(0-10Z3)    =    TIME    ARRAY    (SEC) 

V(0-1023>    =    VOLTS   ARRAY    (VOLTS) 

THE    SUBTRACTION   RESULTS    ARE    THEN   ABLED   TO   BE    PLOTTED   AS    HARD 
COPY    ON    THE    THINK    JET   OR   PEN   PLOTTER.    A    STANDARD   OPTION    IS 
INCORPORATED      FOR   EASE    WHEN   USING    THE    PROGRAM   FREQUENTLY   WITH 
THE    SAME    PARAMETERS.       THE    PROGRAM   ALLOWS    FLEXIBILITY    IN 
CHANGING    THE   PLOT'S   PARAMTERS    BEFORE    THE    HARD   COPY. 

ii2t  luiiiiiirir  iriijz  ti  II  It  itiiarit  iTii  luiirii  2^ 

INITIALIZATION 
iTTT  Tt  uiT  iin  uuTtuiTTrii  iiPniiTTPnPHira      itniitrrtii  n  u  iiiriiitiiTriiiiTijiTTTi  II  null  ii  IT  niriiitiiiiTrir  trmi 

OPTION   BASE    0 

DIM 

DIM 

DIM 

DIM 

DIM 


X_da-ta(  1023)  .Y_da-ta(  1  02  3  )  ,  Y_ax  i5_uni  ts$[  321 
Head_lbl$[G0] , Y_axis_lbl$[ 60] , X_ax is_lbl$[ 32 ] 
X_axis_unit5$[321 . Clear_crt$[ 2 ] ,C$[2] ,Sub_head_lbl$[ 321 
X_char$[8]  ,Vdata$[  15]  ,W$[40]  ,Vda-ta2$[  15]  ,Y_char$[8] 
X 1  (  1  02  3  )  ,  Y  1  (  1 02  3  )  .  XZ  (  1 0Z  3  )  .  YZ  (  1  02  3  ) 

Clear_crt$=CHR$(Z55)&CHR$(75) 

C$=CHR$(Z55)&"K" 


J*     *.      ■.     *.     ■.    ^j  ^^^^^^^^^^_^^L^^    ..      a.      ..      ..      ..      ..      ..     ..      .■     » 

irirTriTTriiiTiTirirTritnaTriiir  u  nniTi 


.    .t«.    .«    .»    J.    »»    «»    »»    ..    .«    ..    ■»    ..    ..    .»    .«    ..    ..    «.    ..    ..    ..    ..    ■»    ».    ■*    ».    ..    .»    Jt 

f  TT  717111  nnirTiTinnn  u  irir  PirnniiTTiTiTniririrTTirTi 

T^TI  TT  TT  ^T^W^W^^^T^T^TT^T^T^ttT^T^T^T^T  FI    II   TT  TI    il   IT   Fl    IT  TI    Fl    11    TI 


GOSUB  Load_file 

INPUT  "DO  YOU  WISH  TO  USE  THE  STANDARD  PARAMETERS^  Y=YES" ,DecZ$ 

IF  DecZ$="Y"  THEN 

GOSUB   S'tand_parapi5 

GOTO   420 
END   IF 

GOSUB  Init ial_parans 
A$="N'" 

GOSUB  Plot_data 

INPUT  "WISH  TO  CHANGE  CURRENT  GRAPHICS  PARAMETERS?  Y=YES".Dec2$ 
IF  Dec2$="Y"  THEN  GOSUB  Change_parans 
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460  INPUT  "ANOTHER  DATR  FILE  PLOTTED  ON  SRME  HRRD  COPY  GRAPH?" ,Decb$ 

470  WHILE  DecG$="Y" 

480  GOSUB  Load_file 

490  INPUT  "LINE  TYPE( 1 =S0LID,4=DfiSH . 5=BIG  DASH , G=DflSH- DOT" .Line.nun 

500  INPUT  "PEN  NUMBER  1  TO  8" .Pen_num 

5t0  GOSUB  Same_graph 

5Z0  GOTO  460 

530  END  UHILE 

540  INPUT  "DO  YOU  WISH  TO  CHANGE  ANY  MORE  PARAMTERS?  Y=YES".DecZ$ 

550  IF  DecZ$="Y"  THEN  GOTO  450 

560  INPUT  "DO  YOU  UISH  TO  DUMP  GRAPHICS?  Y=YES".A$ 

570  IF  A$-"Y"  THEN  GOSUB  Plot_data 

580  INPUT  "DO  YOU  UANT  TO  PLOT  A  DIFFERENT  GRAPH?  Y=YES",Dec7$ 

590  IF  Dec7S="Y"  THEN  GOTO  350 

G00  GOSUB  End_progran 

G10  ! 

630  I               STANDARD  PARAMETERS 

^^      M     /%  I        ^^-^^^^    *-*     A  t     A   a     ^^    A   A       .   A      AA      fc-fc     A   A      A   J       At       *«       ■*      Aa       m    a       a    *      a   a      a  A       Ml      A  A      AM      A   A      J^       .J       **      **      *-       A   A       A  A       *»       A  a      ^      ±|      ^^      **      ■   *      '    •       '-       *-       **       ■-       »■       «■        '*       *'       'A       *-«a.**>       ^^gj^iltjia**       '*       *'       **       *'       ** 

640  I  »4t»#t;itBggtt4Ttt»»«1t«gtTHttitHtntntti^^ 

650  I 

GG0  Stand_params :   ! 

670  Head_lbl$=U$ 

680  SuD_head_lbl$=Vdaia$ 

690  X_char$-="K" 

700  Y_char$="DD.DD" 

710  Y_axis_lbI$="Volts" 

720  X_axi5_lbl$-"Relative  Tine" 

730  Y_axi5_unit5$="Volt5" 

740  X_aA  is_units$  =  "nicrosecond5" 

750  Left=0 

760  Right=19»Tscale 

770  Botton«=-V5cale»4 

780  Top=-Vscale*4 

790  X_tic_5p-Tscale/Z 

800  Y_tic_sp=V5cale/5 

810  X_naj_c-t  =  Z 

8Z0  Y_(na3_ct  =  5 

830  Ma]_tic-Z 

840  Xgndtic-Tscale 

850  Ygndtic-Vscale 

8B0  Xgridna3="1 

870  Ygridmaj-1 

880  Majgrid_s2-Z 

890  Y_incfit  =  Z»V5cale 

900  X_incnt*4»T5cale 

910  Line_nuin*=1 

9Z0  RETURN 

930  I 

HAW  I  iTa HB itnnu nnitnHirnn iiTrna n irnttnir    nniiiriiciTirnn mum iiiriiiiTiiiiiiiTTiiTiiiHiriiHii n nunmr 

950  !              INITIAL  PARAMETERS 

n     f^/&  I          --      --      --      >»      ■-      »*      «»      J-      J«      •*      ■•      »1      »»      »«      «J      ..      ..       .J      *-      -■       •'      J.      ..      .-       ■■      --       ■•        -A      ..      1.      J»      ..       ..       ..      »»      ..       .-       AA       .A       Aa      aA      ..       ..      ^.       ..      ..      ».      A.      ..      ..      ..      JA      »»      IJ      11      ..       .A      AJ       .»      A»      AA      AA     A.A      A-A 

H  hvj  I  HniinnnniipnanPnnnnitniiiiiriiitiinu  tin  it  itirnniiitii  niiTriiiniH  n  uHirjnnrpirn  u  n  TTimniniiiTi 
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I 
ni t ial_parans :    ! 
X_-tic_so=0 
Y_tic_sp=0 
X_naj_ct= 1 
Y_maj_ct= 1 
Maj_t ic=Z 
Xgr  id"l  ic-Q) 
Ygridt ic=0 
Xgr idnaj  =  1 
Ygr  idnaj  =  1 
Majgrid_sz=Z 

INPUT  "GRftPH  HEftDING  LftBEL . 3Z  CHAR.  MftX" .Head_lbl$ 
'GRAPH  SUB  LABEL. 3Z  CHAR.  MAX" , Sub_head_lbl $ 


INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INl^UJT 

Dr-^p 


"X 
"Y 
"X 


970 

980  I 

990 

1000 

1010 

10Z0 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1  100 

1  110 

11Z0 

1  130 

1  140 

1  150 

1  160 

1  170 

1  171 

1  180 

1  190 

1Z00 

1Z10 

1ZZ0 

1Z30 

1Z40 

1Z50 

1ZB0 

1Z70 

1Z80 

1Z90 

1300 

1301 

1310 

13Z0 

1330 

1340 

1350 

1360 

1370 

1 380  Load_f ile; 

1390   ! 

1400 

1410 

14Z0 


"IMAGE  FOR  X  AXIS  NUMERICS. 
"IMAGE  FOP  Y  AXIS  NUMERICS, 
"Y  AXIS  L^^HEL,3Z  CHAR.  MAX" 


AXIS  LABEL 
AXIS  UNITS 
AXIS  UNITS 


8  CHAR. 
8  CHAR. 
Y_ax 1 5_ 
X_ax i5_ 
Y_ax is_ 
X  axis 


MAX"  . 

MAX". 
.IblJ 
lbl$ 
unitsS 
uni-ts$ 


char$ 
char$ 


3Z    CHAR.    MAX" 
3Z    CHAR.    MAX" 
3Z    CHAR.    MAX" 
"UINDOU    LOOKED   AT . LEFT . RIGHT .BOTTOM . TOP" ; 
INPUT    Left  .Right  .Bot-ton, Top 

INPUT    "X    AXIS    TICK    SPACING . DEF  =  0(  NO    TI CKS ) " . X_t ic_sp 
"Y    AXIS    TICK    SPACING,DEF=0(NO   TICKS  )". Y_t ic_sp 
"NUfI    TIC    INTVLS    BTUN   MAJ   TICS   X-DIR  ,DEF=  1"  ,  X_naj_ct 
"NUM    TIC    INTVL5   BTUN   MAJ    TICS   Y-DIR,DEF= 1" . Y_naj_ct 
"SIZE    OF    MAJOR    TICKS .DEF-Z" .Maj_t ic 
"GRID   X-TIC    SPACING.DEF=0<NO    TICKS  )". Xgridt ic 
INPUT    "GRID   Y-TIC    SPACING  .DEF=0(  NO    TICKS  )".  Ygndt  ic 
INPUT    "NUM    TICKS   BTUN   MAJOR   TICKS    Y-DIR .DEF= I" . Ygr idnaj 
"NUM    TICKS    BTUN   MAJOR   TICKS    X-DIR.DEF=  1"  .  Xgndnaj 
"SIZE    OF    MAJOR  GRID   TICKS ,DEF=Z" .Majgrid_sz 
"Y-AXIS    INCREMENT" .Y_incnt 
"X-AXIS    INCREMENT". X_incmt 
LINE    TYPE( 1=S0LID,4=DASH,5-BIG   DASH,6=DASH-O0T ) " ; 


INPUT 
INPUT 
INPUT 
INPUT 
INPUT 


INPUT 

INPUT 

INPUT 

INPUT 

DISP 

INPUT   Line_nun 

! 

RETURN 

n  n  "  "  "  "  n  Tib  niliriInTIHII"""~II""  "  ■ 

!  tOAD  DATA  FILE 


I 


INPUT  "INPUT  NAME  OF  DATA  FILE  l".Vdata$ 
INPUT  "INPUT  NAME  OF  DATA  FILE  Z  ".VdataZJ 
ASSIGN  eOfl  TO  VdataS 
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1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
16Z0 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
17:0 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 


ENTER  @Df 1 ;U$,Tscale,V5cale,X1(  »  > .Y1(  •  ) 

ASSIGN  @Df1  TO  » 

ASSIGN  @DfZ  TO  VdataZ$ 

ENTER  @DfZ;U$,Tscale,V5cale.XZ(»).YZ(» ) 

ASSIGN  @Df2  TO  • 

MAT  K_data=  XI 

MAT  Y_data=  Y1-Y2 

PRINT  "DATA  FILES  READ" 

BEEP 

RETURN 


«ii    tA   ii   *■    **    **    **    **    **    ■*    **    ■*    *■    **    ** 
itiiiiiiiriip  n  niTBini  ii  ii 


*■     **     ■*     *■     ■*     **     **     ■*      **     **     *  '      *fc     ■*      ■■     **     fci     ii      ii      ■'      '  *     *■     **     J*      *-      M.  ±     **      *■      >*     *■     M.  i     AK     fc*      fcj     Aa     **     t 

iiiiitiiiiiriinii  n  imiiiriiirniin  irii  n  n  TTirit  it  ititit  ii  it  ititiii 

PLOT   DATA 

gHnitu  n nuitirirTTlI  ullltltitltultitirua u Tlitltp  irnTi  uTTTrilirilllltlllt 

^T  TT  TT  IT   f  f    9t  ^T  ^T^T^T^W^TT^T^T^T^T^T^T^^^T^T^T^^  ^T  ^T^T  TT  ^T^T  ^T^^  TT  ^TTT^T^^^T  PT   If    II    fl    II    11    ff  T^^T^T^t^T^tTF^W^  TT  TT^T^^TT  ^T^TTTT  TT   TT  TT  TT 


Plot_data:   ! 
I 

IF  A$="Y"  THEN 

OUTPUT  KBO;Clear_crt$; 

GINIT 

PLOTTER  IS  CRT, "INTERNAL" 

GOTO  1770 
END  IF 

INPUT  "IS  THIS  TO  BE  A  HARDCOPY  PLOT?  Y=YES".Decl$ 
IF  Df'cl$  =  "Y"  THEN 

GINIT 

PLOTTER  IS  705."HP6L" 

GOTO  1770 
ELSE 

OUTPUT  KBD;Clear_crt$; 

GINIT 

PLOTTER  IS  CRT. "INTERNAL" 

GOTO  1770 
END  IF 
GRAPHICS  ON 
PEN  1 

X_gdu_nax=100»MAX( 1 .RATIO) 
Y_gdu_nax-100*MAX( 1 , 1/RATIO) 
LORG  5 
CSI2E  4 

MOVE  .48»X_gdu_max , .93»Y_gdu_nax 
LABEL  Head_lbi$ 

MOVE  .48»X_gdu_nax  ,  .89*Y_gdu_iiax 
CSIZE  3 

LABEL    Sub_head_lbl$ 
DEG 

LDIR   90 
CSIZE    4 
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1910 

19Z0 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

Z000 

2010 

Z0Z0 

Z030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

21  10 

2120 

ZI30 

2140 

2  150 

216C 

2  170 

2  180 

2190 

2200 

ZZ10 

2220 

2230 

2240 

2250 

2260 

2270 

Z280 

2290 

2  300 

2310 

2320 

2  330 

2  340 

2  350 

2  360 

2  370 

2  380 


LRBEL 
MOVE 
CSI2E 
LftBEL 


15*Y  Qdu  max 


l9»Y_gdu_nax 


MOVE  .04»X_gdu_nax , .55»Y_gdu_nax 

LABEL  Y_axi5_lbi$ 

LORG  5 

LDIR  0 

MOVE  .48*X_gdu_nax 

CSIZE  4 

X_axis_lbl$ 
48*X_gdu_max 
3 

"(•';X_axis_units$;"  )" 
MOVE  .70»X_gdu_max , .74»Y_adu_nax 

VirUPORT  .  15»X_gdu_nax.  .  85»X_gdu_nax  , .  25»  Y_gdu_nax  , .  86»Y_gdu_fiax 
UINDOU  Left .Right .Bottom. Top 

ftXES  X_t ic_sp , Y_t ic_5p .Left, Bottom . X_maj_ct , Y_maj_ct . Maj  _t  ic 
AXES  X_t  ic_sp , Y_t  ic_5p .Right .Top ,X_maj_ct , Y_maj_ct ,Ma3_t ic 
IGRID  Xgridt  ic,  Ygridt  ic  .Left  .Bottom. Xgndmaj  .Ygndmaj  ,Majgrid_52 
CLIP  OFF 
CSIZE  3, .75 
LORG  6 
FOR  I=Left  TO  Right  STEP  X_incmt ! LABELING  OF  X  AXIS 

MOVE  I .(Bottom-. 04»( Top) ) 

LABEL  USING  X_charJ; I / 1 . E-G 
NEXT  I 
LORG  8 
FOR  I=Bottom  TO  Top  STEP 

MOVE  (Left-.4»T5cal.'j)  . 

LABEL  USING  Y_char*;I 
NEXT  I 
PENUP 
PEN  2 

LINE  TYPE  Line_num 
CLIP  ON 
FOR  1=0  TO  1023!PLOTTING  OF  DATA 

PLOT  X_data( I ),Y_data( I) 
NEXT  I 
PENUP 
CLIP  OFF 
IF  A$="Y' 
RETURN 


Y_incmt! LABELING  OF  Y  AXI: 
I 


THEN  DUMP  GRAPHICS  «701 


1t1t1t1t1t  tltttttltttttltHttllltttttttlHtttttttltlt 
SAME    GRAPH 


I 


Itttttttltttttt 


W  W  11  TT  W  HIT 


Same  graph: 
PEN  Pen_num 
LINE  TYPE  Line. 
CLIP  ON 


num 
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I 


Z390 

Z400 

Z410 

2420 

Z430 

2440 

2450 

Z460 

2470 

2480 

2490 

Z500 

2510 

2520 

2530 

2540 

2550 

2  5G0 

2570 

Z580 

2590 

2600 

2610 

Z6Z0 

ZB30 

2640 

2E50 

Z6G0 

2670 

2680 

2630 

2700 

2710 

Z720 

Z730 

Z740 

Z750 

Z760 

Z770 

Z780 

2750 

Z8«0 

28!0 

Z8Z0 

Z830 

Z840 

Z850 

Z860 


FOR  1=0  ro  !0Z3 

PLOT  X_data(I  ),Y_data(  I  ) 
NEXT  I 
CLIP  OFF 
PENUP 
RETURN 

CHANGE  GRAPHICS  PARRMETERS 

aj      .1      m  M     *i      .*      >*     KM      *■     Aj    ^^    *■      ■*      fcfc      fc*     *■      fcj      **      tj      aj     at      '*      litj     t  t      kt      ii     **      **      **     "'      *■      tj      1  t     *  ■     aa     **     ij      *j      afc     *>     **     *■      *■      **     *j     **     **     *»      **      ■*     ■*      >*      ****     *■*■     A*     ***j     *AA*jfc      La     aa. 

all  naiiiiniriiiinTritirirnnirnnTTgTTnTi  nn  n  niiaiiHitnirattnitiittnnii  n  itgiinirH  nn  mumniiH  it  im 

I 

Change_paran5:   • 
OUTPUT  KBD;Clear_crt$; 
PRINT  "SELECT  FROM  THE 


FOLLOUING  TO  CHANGE- 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
P^INT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
I 

INPUT 
IF 


"CASE  1 :  GRAPH  HEADING" 

"CASE  Z:  Y  AXIS  LABEL" 

"CASE  3:  X  AXIS  LABEL" 

"CASE  4:  Y  AXIS  UNITS" 

"CASE  5:  X  AXIS  UNITS" 

"CASE  6:  UINDOU  UNITS" 

"CASE  7:  X  AXIS  TICK  SPACING" 

"CASE  8:  Y  AXIS  TICK  SPACING" 

"CASE  9:  TICK  INTVL  STUN  MAJOR  TICKS  IN  X  DIRECTION" 

"CASE  10:  TICK  INTVL  BTUN  MAJOR  TICKS  IN  Y  DIRECTION" 

"CASE  11:  SI7E  OF  MAJOR  TICKS" 

"CASE  IZ:  X  TICK  GRID  SPACING" 

"CASE  13:  Y  TICK  GRID  SPACING" 

"CASE  14:  GRID  TICK  INTVL  BTUN  MAJOR  TICKS 

"CASE  15:  GRID  TICK  INTVL  BTUN  MAJOR  TICKS 

"CASE  16:  MAJOR  GRID  TICK  SIZE" 

"CASE  17:  Y  INCREMENT" 

"CASE  18:  X  INCREMENT" 

"CASE  19:  LINE  TYPE" 

"CASE  20:  GRAPH  HEADING  SUB  LABEL" 

"CASE  Zl:  X  AXIS  NUMERIC  IMAGE" 

"CASE  Z2:  Y  AXIS  NUMERIC  IMAGE" 


IN 
IN 


DIRECTION" 
DIRECTION" 


"DO  YOU  DESIRE  TO  CHANGE  ANY  PARAMETERS?  Y-YES".D«c4$ 
Dec4$-"Y"  THEN 
INPUT  "ENTER  PARAMETER  NUMBER" , Change.num 
ELSE 

GOTO  3560 
END  IF 

SELECT  Change_nun 
CASE  1 

PRINT  "CURRENT  HEADING  LABEL" .Head_lbl$ 
INPUT  "NEW  HEADING  LABEL" .Head_IbI$ 
CASE  Z 
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Z870        PRINT  "CURRENT  Y  ftXIS  LRBEL" . Y_ax is_lbl$ 

Z880        INPUT  "NEU  Y  ftXIS  LRBEL" , Y_ax is_lbl $ 

Z890  CRSE  3 

2900        PRINT  "CURRENT  X  RXIS  LABEL" , X_ax is_lbl$ 

Z910        INPUT  "NEU  X  AXIS  LABEL" . X_ax is_l bl $ 

Z9Z0  CASE  4 

2930        PRINT  "CURRENT  Y  AXIS  UNITS" . Y_ax is.unitsS 

2940        INPUT  "NEU  Y  AXIS  UNITS" . Y_ax is_units$ 

Z950  CASE  5 

Z9B0        PRINT  "CURRENT  X  AXIS  UNITS" . X_ax is_units$ 

2970        INPUT  "NEU  X  AXIS  UNITS" . X_ax is_units$ 

Z980  CASE  B 

Z990        PRINT  "CURRENT  UINDOU"  .Left  .  Right  ,Bot-t  on  ,  Top 

3000        INPUT  "NEU  UINDOU"  .Left  .Right  .Bottopi. Top 

3010  CASE  7 

30Z0        PRINT  "CURRENT  X  AXIS  TIC  SPACING" . X_t ic_sp 

303(d        INPUT  "NEU  X  AXIS  TIC  SPACING"  .  X_t  ic_5p 

3040  CASE  8 

3050        PRINT  "CURRENT  Y  AXIS  TICK  SPACING" . Y_t ic_5p 

30B0        INPUT  "NEU  Y  AXIS  TIC  SPACING" . Y_t ic_5p 

3070  CASE  9 

3080        PRINT  "CURRENT  TIC  INTVL  BTUN  MAJ  TICS  X  DIR" . X_naj _ct 

3090        INPUT  "NEU  TICK  INTVL  BTUN  MAJOR  TICKS  X  DIR" . X_naj _ct 

3100  CASE  10 

3110        PRINT  "CURRENT  TICK  INTVL  BTUN  MAJOR  TICKS  Y  DIR" . Y_maj_ct 

31Z0        INPUT  "NEU  TICK  INTVL  BTUN  MAJOR  TICKS  Y  DIR" . Y_maj_ct 

3130  CASE  n 

3140        PRINT  "CURRENT  MAJOR  TICK  SIZE" , Maj_t ic 

3150        INPUT  "N£U  MAJOR  TICK  SIZE" . Maj_t ic 

31B0  CASE  12 

3170        PRINT  "CURRENT  X  TICK  GRID  SPACING" .X_grid_t ic 

3180        INPUT  "NEU  X  TICK  GRID  SPACING" .X_grid_tic 

3190  CASE  13 

3Z00        PRINT  "CURRENT  Y  TICK  GRID  SPACING" .Y_grid_t ic 

3Z10        INPUT  "NEU  Y  TICK  GRID  SPACING" . Y_gr id_t ic 

3Z20  CASE  !4 

3Z30        PRINT  "CURRENT  GRID  TIC  INTVL  BTUN  MAJ  TICS  X  DIR"  .X_grid_Piaj 

3240        INPUT  "NEU  GRID  TIC  INTVL  BTUN  MAJ  TIC  X  OIR"  ,X_Qrid_inaj 

3250  CASE  15 

32G0        PRINT  "CURRENT  GRID  TIC  INTVL  BTUN  MAJ  TICS  Y  DIR"  .  Y_grid_r«aj 

3Z70        INPUT  "NEU  GRID  TIC  INTVL  BTUN  MAJ  TIC  Y  DIR" , Y_grid_naj 

3Z80  CASE  IB 

3290  PRINT    "CURRENT   MAJOR   TICK    SIZE" .Ma]_grid_sz 

3300  INPUT    "NEU   MAJOR   TICK    SIZE" . Maj_grid_52 

3310  CASE    17 

33Z0  PRINT    "CURRENT    Y    INCREMENT"  ,  Y_incr>it 

3330  INPUT    "NEU   Y    INCREMENT"  .  Y_incPit 

3340  CASE    18 


I 


3350 

33G0 

3370 

3380 

3330 

3391 

3400 

3410 

34Z0 

3430 

3440 

3450 

34G0 

3470 

3480 

3490 

3500 

3510 

35ZO 

3530 

3540 

3550 

35B0 

3570 

3560 

3530 

3600 

3510 

3G20 

3G30 

3B40 

3650 


PRINT  "CURRENT  X  INCREMENT" . X_incmt 
PRINT  "NEU  X  INCREMENT" ,X_incmt 
CASE  19 

PRINT  "CURRENT  LINE  TYPE" . Line_nun 

DISP  "NEU  LINE  TYPE(  1  =S0LID.4=DflSH  ,5=BIG  DflSH.G=DflSH-DOT"  j 
INPUT  Line_nun 
CASE  Z0 

PRINT  "CURRENT  HEADING  SUB  LABEL" . Sub_head_lblS 
INPUT  "NEU  HEADING  SUB  LABEL" , Sub_head_l bl$ 
CASE  Zl 

PRINT  "CURRENT  X  AXIS  NUMERIC  IMAGE" . X_char$ 
INPUT  "NEU  X  AXIS  NUMERIC  IMAGE" . X_char$ 
CASE  ZZ 

PRINT  "CURRENT  Y  AXIS  NUMERIC  IMAGE" . Y_char$ 
INPUT  "NEU  Y  AXIS  NUMERIC  IMAGE" . Y_char$ 
END  SELECT 

INPUT  "DO  YOU  UISH  TO  CHANGE  ANOTHER  PARAMETER?  Y=YES" ,Dec5$ 
IF  Dec5$="Y"  THEN 

GOTO  Z510 
ELSE 

GOSUB  Plot_data 
END  IF 
RETURN 


IlTlTtTrH  niTH  aalTB  uauQlIItlru  XtiltZTllsT^XittaiitiutiuuTXiiitu  o  a 

END    THL"    PROGRAM 

B  nniTH  lIItllTTnTra  H  KH  ITTTnllllTtB  allTrilirTr  nirinrilll  a  a  a  UmraTilinTI  il  IT  n  murTTTrB  tITTTi  g  a 


End_progran:   ! 

GRAPHICS  OFF 

PRINT  "END  OF  PROGRAM" 

END 
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